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We investigate phase structure of the D{> 4)-dimensional gauged Nambu-Jona-Lasinio (NJL) 
model with 5{= D — 4) extra dimensions compactified on TeV scale, based on the improved ladder 
Schwinger-Dyson (SD) equation in the bulk. We assume that the bulk (dimensionless) running 
gauge coupling in the SD equation for the SU{Nc) gauge theory with Nf massless flavors is given 
by the truncated Kaluza-Klein effective theory and hence has a nontrivial ultraviolet fixed point 
(UVFP), resulting in the walking coupling. We find the critical line in the parameter space of 
two couplings, the gauge coupling (value fixed at the UVFP) and the (dimensionless) four-fermion 
coupling, which is similar to that of the gauged NJL model with fixed (walking) gauge coupling 
in four dimensions. It is shown that in the presence of such walking gauge interactions the four- 
fermion interactions become "nontrivial" and "renormalizable" even in higher dimensions, similarly 
to the four- dimensional gauged NJL model. Such a renormalizability/nontriviality holds only in 
the restricted region of the critical line ( "nontrivial window" ) with the gauge coupling larger than 
a non-vanishing value ("marginal triviality (MT)" point), in contrast to the four- dimensional case 
where such a renormalizability holds for all regions of the critical line except for the pure NJL point 
(without gauge coupling). In the nontrivial window the renormalized effective potential yields a 
nontrivial interaction which is conformal invariant. The existence of the nontrivial window implies 
"cutoff insensitivity" of the physics prediction in spite of the ultraviolet dominance of the dynamics. 
In the formal limit D — + 4, the nontrivial window shrinks to the pure NJL point but with a 
nontrivial condition which coincides with the known condition of the renormalizability/nontriviality 



of the four-dimensional gauged NJL model. 
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PACS numbers: 11.15.Ex,11.10.Kk,11.25.Mj,12.60.Rc 



I. INTRODUCTION 



Growing attention has recently been paid to the model 
buildings based on the scenario of extra dimensions com- 
pactified with TeV scale [l|, |3]- Along such a scenario 
the top quark condensate (Top-Mode Standard Model: 
TMSM) l a B, la lo l was reconsidered by several au- 
thors 0, Hj S Il3i US nil in the versions with compact- 
ified extra dimensions, where the ad hoc four-fermion 
couplings in the original TMSM may be replaced by 
the strong attractive forces due to Kaluza-Klein (KK) 
modes of the bulk Standard Model (SM) gauge bosons. 
In particular, Arkani-Hamed, Cheng, Dobrescu and Hall 
(ACDH) proposed a version of the TMSM where the 
third generation quarks/leptons as well as the SM gauge 
bosons are put in the bulk, while other generations in the 
brane. In a series of papers [T^, UM ^^ investigated the 
full gauge dynamics {without four-fermion interactions) 
in the bulk, based on the ladder Schwinger-Dyson (SD) 
equation. We further studied lz\ phenomenological im- 
plications on the ACDH scenario and found that D = 8 
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( four extra dimensions) will be a viable case which pre- 
dicts a correct top mass, nit = 172 — 175 GeV, and a 
characteristic Higgs mass, niH = 176 — 188 GeV, from 
the requirement that only the top can condense ( "Topped 
MAC" or tMAC requirement). 

Actually, the most prominent feature of the gauge the- 
ories with extra dimensions is that the theories become 
strongly coupled, with the coupling growing by the power 
running in the ultraviolet region beyond the compactifi- 
cation scale. Let us take a class of SU{Nc) gauge the- 
ories with massless Nf flavors. It was pointed out Jm\ 
(see also Il3| |l3, 0) that the bulk gauge coupling of 
such gauge theories can have a (nontrivial) ultraviolet 
fixed point (UVFP) within the truncated KK effective 
theory :2j . Although existence of such a UVFP is some- 
what subtle in the lattice studies [ifil llTl Il8l | and other 
nonperturbative approach jl9|, unless it is an artifact of 
the truncated KK effective theory, it implies existence 
of a renormalizable interacting gauge theory even in di- 
mensions higher than four. It is in sharp contrast to 
the conventional view (based on the perturbation) that 
higher dimensional gauge theories are nonrenormalizable 
and trivial. Once we assume that such a nontrivial UVFP 
does exist beyond the truncated KK effective theory, the 
theory behaves near the UVFP as a walking gauge the- 
ory 20, 21, 22, 23, 24 1 ha ving a large anomalous dimen- 
sion 7m = D/2 ~ 1 jiri HJ. In the four-dimensional 
case, as was pointed out by Bardeen et al. |23, the four- 



fermion interaction becomes a marginal operator due to 
the large anomalous dimension, 7j„ ~ 1, of the walking 
gauge theory, and hence can be mixed with the gauge in- 
teraction. Then the phase structure of the walking gauge 
theories should be analyzed in a larger parameter space in- 
cluding the four-fermion interactions as well as the gauge 
interactions. 

The gauged Nambu-Jona-Lasinio (NJL) model, a 
model of gauge theory plus four-fermion interaction, was 
first studied by Bardeen, Leung and Love !23| within 
the ladder SD equation for the fixed gauge coupling in 
four dimensions. The critical line of the four-dimensional 
model for the fixed gauge coupling was discovered by 
Kondo, Mino and Yamawaki |2a | and independentl y by 
Appelquist, Soldate, Takeuchi and Wijewardhana |27l |. 
Since then the model has been offering many interesting 
applications for physics of dynamical symmetry breaking, 
such as the walking technicolor 113, iHlHlslill , the 
stron g ET C technicolor |2 qj 2 3 /topcolor-assisted techni- 
color^, and the TMSmHMJ1||] or topcolor jp], etc. 
(See for reviews Refs. |23,l23,l2J|)- 

The characteristic dynamical feature of the four- 
dimensional gauged NJL model is the large anomalous 
dimension |23|, 1 < 7™ < 2, for the fixed (walking) 
gauge coupling, and the finiteness of all (renormalized) 
couplings at composite level as well as induced Yukawa 
coupling of the fermion to the composite bosons in the 
continuum limit (infinite cutoff limit). This implies the 
(nonperturbative) "renormalizability" or "nontriviality" 
of the model, as was pointed out in Ref . |3a, [Sg for the 
("moderately") walking gauge coupling, with A > 1 ^, 
and was subsequently shown more systematically for the 
fixed (walking) gauge coupling {A ^^ 00) in Ref. 13^ a s 
well as for the moderately walking one in Ref. |37l 133 . l39l| . 
For the fixed gauge coupling, such a phenomenon takes 
place along the whole critical line except for the pure NJL 
point (without gauge coupling and "7^ = 2" ) where the 
above couplings vanish logarithmically in the continuum 
limit and hence the theory becomes nonrenormalizable 
or trivial. 

Inclusion of four-fermion interactions into the gauge 
theories has not been done so far in the context of the 
extra dimension scenario ^. The analysis of the phase 
structure of the gauged NJL model with extra dimensions 
for such walking gauge theories may be useful not only 
for the TMSM but for other strong coupling theories such 
as the bulk technicolor, various composite models, etc.. 



1 A = QC'i/i-b), where Cp = (N^ - l)/2Nc and -b = {UNc - 
2Nf)/3. The parameter A measures the speed of running of the 
running gauge couphng 12% h3j| : A < 1 is the normal running as 
in the QCD, while yl S> 1 is the walking coupling and A — > 00 
corresponds to the fixed ( "standing" limit of the walking) gauge 
coupling. Combined with the asymptotic freedom, A > 1 reads 
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■ Nc < f ATc 



In this paper we investigate phase structure of the 
gauged NJL model in the D{> 4)-dimensional bulk, 
based on the improved ladder SD equation for the bulk 
fermion propagator, assuming that the (dimensionless) 
bulk running gauge coupling is given by the truncated 
KK effective theory and hence almost fixed at the non- 
trivial UVFP in a wide energy range above the compact- 
ified energy scale. As far as the cutoff energy A is large, 
the effect of the compactified scale is negligible and the 
SD equation is well approximated by that with the fixed 
coupling [13, lU ■ Thus we take the gauge coupling be 
fixed at the value of UVFP for simplicity. Here we should 
stress that our results in this paper crucially depend on 
the existence of such a UVFP. 

By solving the SD equation we find the critical line in 
the two-dimensional parameter space {no^g) of the cou- 
pling strength kd (essentially the squared gauge coupling 
fixed at the UVFP value) and the (dimensionless) four- 
fermion coupling g: 
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crit _ 2 



kd/kS" 



(LI) 



which exists only for < kd < ^v^J'*, where k™* is the 
critical coupling strength of the pure gauge theory with- 
out four-fermion interaction |l(l lll| . The critical line 
separating the spontaneously broken phase (with nonzero 
dynamical mass M{^ 0) of the fermion) and the unbro- 
ken one (with zero dynamical mass) resembles that in 
the four-dimensional gauged NJL model with fixed gauge 
couphng [2^ [23 except for the prefactor D/2— 1 instead 
of 1. Note, however, that we actually consider in this 
paper a class of SU{N^) gauge theories (with Nf fiavors) 
which are asymptotically free having a trivial UVFP in 
four dimensions, kd — > (D — *■ 4), but do have the non- 
trivial one in higher dimensions, kd 7^ (D > 4) whose 
value is uniquely determined once the model parameters 
such as Nc and Nf are specified. Each point on the 
above critical line thus corresponds to a single gauged 
NJL model with a different gauge theory with different 
iVc, Nf. We cannot freely adjust its value. 

We further discuss the phase structure of the model, 
following the manner in the four-dimensional case |37l| . 
We find that the decay constant F^^ of the composite 
Nambu-Goldstone (NG) boson tt becomes finite in the 
continuum limit [A/M -^ 00), if Nc and Nf are arranged 
so as to satisfy the condition 



Kd < Kd < K^' , (nontrivial window). 



where k^ (7^ 0) is given by 



,MT 



1 



1 



{D/2-iy 



(1.2) 



(1.3) 



^ The pure NJL mode l w i th extra dimensio ns h as been already 
studied, e.g. in Refs. E3.li3l- See also Ref. 0. 



We may call k^"^ as a "marginal triviality (MT)" point, 
where the decay constant F^ of the composite NG bo- 
son TT diverges logarithmically, F^ ~ In A, in the con- 
tinuum limit A/Af — s- cx) and hence resembles the pure 



NJL point (zero gauge coupling) in the four-dimensional 
case. Then in the nontrivial window we find a sur- 
prising result that even in higher dimensions the four- 
fermion theory in the presence of gauge interactions be- 
comes "renormalizable" and "nontrivial" in the sense 
that the induced (renormalized) Yukawa coupling of the 
fermion to the NG boson and the (renormalized) cou- 
plings among tt and composite scalar a all remain finite 
in the continuum limit (A/Af -^ cx)). This is analogous 
to the renormalizability and nontriviality of the four- 
dimensional case (35, 36, 3'3,|33,|33 except that the MT 
point in four-dimensional case corresponds to the pure 
NJL point, K^'^ ^ for £) ^ 4. We demonstrate these 
facts by explicit computation of F^, the Yukawa vertex 
Fg, the scalar propagator Da-, and the effective potential 
y(a,7r). 

Here we should emphasize one distinct feature of the 
inclusion of the four-fermion interaction into the gauge 
theory in the bulk with compactificd extra dimensions: 
In the case of pure gauge theories (without four-fermion 
interactions), the "nontrivial" theory defined at the crit- 
ical point K^'* is only formal, since the gauge coupling 
kd at the UVFP is not a continuous parameter but a dis- 
crete one depending on Nc and Nf and hence cannot be 
fine-tuned arbitrarily close to the critical value [ifl [ll| . 
It means that even if we assume existence of the UVFP, 
the pure gauge theory in higher dimensions can only be 
renormalizable in a formal sense. Thus there is no finite 
theory for pure gauge theories. On the other hand, in 
the gauged NJL model we do have a continuous parame- 
ter, the four-fermion coupling g, and hence the nontrivial 
theory can be defined by fine-tuning g arbitrarily close to 
the criticality of the four-fermion coupling at the critical 
line Eq. (|I.ip . Note that for kd > k™* where no critical- 
ity of the four-fermion coupling exists, there is no finite 
theory even in the gauged NJL model. 

For the theories in the nontrivial window, k^"^ < kd < 
K^'' , we perform explicitly the renormalization of the ef- 
fective potential in a way similar to the four-dimensional 
case given in Ref.l37|. The renormalized four-fermion 
couphng gn also has a UVFP gR{oo) — g'^''*. We then 
find that the theory has a large anomalous dimension at 
the UVFP gfl(oo): 



Ir, 



^'^' 



1 + J1-kd/h' 



crit 
D 



(1.4) 



which takes the same form as that of the four-dimensional 
model 28, !37] when D in the prefactor is set to be Z? = 4. 
Since D/2 > j,n > D/2 — 1, the dynamical dimension of 
the four-fermion interaction, dim(^?/')^ = 2[(Z3 — 1) — 7,„], 
reads D - 2 < dim(f/'V')^ < D for k^'^ < kd < Kg'*, 
and such an operator is a relevant operator similarly 
to the four-dimensional gauged NJL model |37]. Thus, 
as we were motivated by the four-dimensional case, we 
conclude that the phase structure of the gauge theory 
with extra dimensions which behaves as a walking the- 
ory should also be studied with inclusion of the four- 
fermion interaction. As expected, in the renormalized ef- 



fective potential there exist interactions of an and tt^ for 
K^""" < Kjj < K™* (nontriviality), which are conformal 
invariant. The renormalizability/nontriviality implies a 
remarkable fact that the physics prediction is cutoff in- 
sensitive, notwithstanding the ultraviolet dominance of 
the dynamics of both the walking gauge and the four- 
fermion interactions strongly coupled in the ultraviolet 
region, which may be useful for model buildings. 

We further study RG flows of each theory in the non- 
trivial window governed by a UVFP (kd, 5fl(oo)) cor- 
responding to a set of {Nc,Nf) (each point lying on the 
critical line g"^^). Then we find that our renormalization 
is consistently performed even off the UVFP. 

For < Kjj < K^"^ , on the other hand, we find 
that the gauged NJL model becomes trivial and non- 
renormalizable in spite of the fact that we can fine-tune 
the four-fermion coupling to make finite the dynami- 
cal mass of the fermion. The "renormalized effective 
potential" actually has no interaction terms of an and 
TTj^. This renormalization would lead to the "anoma- 
lous dimension" 7^ = -j and dim(V'V')^ = _D — 2 for 



< Kd < K 



MT 



However, our renormalization is incom- 



plete for < Kd < ^d"'"- There exists another marginal 
operator 9m('0'0)^^^(V'V') whose dynamical dimension is 
given by dim(9M('0^)9*^(^/''0)) = D (marginal opera- 
tor!) in our renormalization scheme. 

We further discuss in Sec. lVIII the formal limit D ^ 4, 
in which k^"^ -^ and kd — > and hence the UVFP's for 
both regions of the nontrivial window k^"^ < kd < k^'* 
and trivial region < kd < k^"^ shrink to a single point 
of the pure NJL point. Nonetheless, the nontrivial win- 
dow for Z? ^ 4 coincides with the condition of the renor- 
malizability/nontriviality of the four-dimensional gauged 
NJL model with (moder ately ) walki ng g au g e coupling 
characterized by A > 1 [3^ ISo. 1371 133 . l39| . which is 
clearly distinguished from A < 1 corresponding to the 
trivial region < kd < k^"^ . Thus the nontrivial win- 
dow is not a peculiarity of the extra dimensions but ac- 
tually exists already in the four-dimensional gauged NJL 
model in a rather sophisticated manner. 

The paper is organized as follows: In Sec. ^we give 
our model setting, Lagrangian and the running gauge 
coupling constant within the truncated KK effective the- 
ory. In Sec. IIHI we study phase structure of the model, 
based on the simplest approximation, the bifurcation 
technique |4J|, to the SD equation. We find the criti- 
cal line, and calculate simplified scaling relation, scalar 
boson propagator and the decay constant F^^ of the NG 
boson TT. We also find the MT point where F^ diverges 
logarithmically in the continuum limit. In Sec. II VI we in- 
corporate subleading terms in the asymptotic region by 
more involved approximations than the bifurcation, "lin- 
earized a ppr oximation" 32, 42J and "power expansion 
method" |43l | . The sub-dominant term is required for the 
effective potential. In Sec.|3we calculate the effective po- 
tential. In Sec. lVII we discuss the renormalization making 
finite the induced Yukawa coupling, the effective poten- 
tial and the scalar boson propagator. We give detailed 



discussions on the RG flow in the nontrivial window. 
Sec. IV ill is devoted to Summary and Discussion: In par- 
ticular we discuss the D -^ 4 limit. In Appendix IM we 
give another detailed calculations of the composite scalar 
propagator. AppendixiBlandOare for details of the "lin- 
earized approximation" and "power expansion method" , 
respectively. Appendix IDI is for more discussions on the 
nonrenormalizability of the region < kd < i^^"^ ■ 



II. THE MODEL 

Let us start with the renormalization group properties 
of the gauge interaction in the D — 4 + S dimensional 
bulk. Extra S ~ [D ~ A) spatial dimensions are assumed 
to be compactified at scale R~^. The negative mass di- 
mension of the bulk gauge coupling implies strong inter- 
action at high energy. Thus the higher loop effects some- 
how need to be taken into account. The task, however, 
immediately meets serious trouble with the nonrenormal- 
izability of the bulk gauge theory D > A. We therefore 
define the bulk gauge theory by using the "truncated 
Kaluza-Klein (KK)" effective theory 0], which allows us 
to calculate the loop effects within 4D effective field the- 
ory. A remarkable feature of the truncated KK effective 
theory is the existence of the non-trivial ultraviolet fixed 
point (UVFP) at least for its MS couphng ^10, J^Ollli. 
We briefly review the properties of the renormalization 
group equation (RGE) in this model. 

After the decomposition of the bulk gauge field into 
its KK-modes, the running of the four-dimensional gauge 
couphng 34£) can be described by the RGE, 



and thus 



d^ 



(II.l) 



where we assumed the renormalization scale /i is suffi- 
ciently larger than the compactification scale i?^^, /^ ^ 
R~^. The RGE coefficient b' is given by 

6'^-^iV, + ^2^/2Ar,, (II.2) 

for SU{Nc) gauge theory with Nf bulk fermions. Nkk{iJ') 
stands for the number of KK modes below the renormal- 
ization scale fj,. We use an approximation 



1 ^*/' 
^^' 2"r(l-h(5/2) 



(II.3) 



The factor 1/2" in Eq. (|II.3p arises from the orbifold com- 
pactification on TV^s"- (See Refs. mill.) The four- 
dimensional gauge coupling g^j^ can be matched with the 
dimensionful bulk gauge coupling gu_^g\jj, 

gl^^,)o = (2^i?)^g2^/2". 

It is convenient to define the dimensionless bulk gauge 
couphng g, 

9 = 5(4+(5)£)A* I 



,2 {2TTRfiy 2 f ^ 

? W = — t;;, — 94dW- 



(11.4) 



Combining Eqs. (|II.ip . I|II.3|I . and (|II.4|I . we obtain RGE 
for g 10] 



V^2'- 



1 



-]n^uAb'g\ (A*»i?-^), (II.5) 



with r^NDA being the loop factor of naive dimensional 
analysis (NDA) in D dimensions. 



n 



1 



NDA 



(47r)^/2r(Z3/2) ■ 
The RGE JUS) leads to a UVFP g* [13, 

1 



5*^NDA 



-{l + 2/S)b'' 



(II.6) 



(II.7) 



for b' < 0. 

So far we have shown the UVFP only within the 
truncated KK effective theory at one-loop level. Does 
such a UVFP really exist beyond the approximation we 
adopted? This question is, of course, a highly non- 
perturbative problem 16, 17, 18.j 19] and extremely dif- 
ficult to be answered. Instead of solving this difficult 
problem, in this paper, we simply assume the existence 
of the UVFP and address a hopefully easier but still quite 
exciting question: How does the bulk gauge theory be- 
have with such a UVFP? 

Such a theory should possess an approximate confor- 
mal invariance and was shown to have a large anoma- 
lous dimension jm = D/2 — 1 [ifl uM- The situation 
has a strong resemblance to the four-dimensional walking 
gauge theories |2y, |2l|, |23, |2J] , in which the gauge cou- 
pling is assumed to be on the nontrivial UVFP and the 
fermion bilinear operator tprp acquires a large anomalous 
dimension 7^ = 1. It is known that the Nambu-Jona- 
Lasinio (NJL) type four-fermion interaction becomes rel- 
evant in such a walking gauge dynamics. The walking 
gauge theory was then analyzed in an extended coupling 
space including the NJL type four-fermion interaction 
(gauged NJL model) '25', ^ lUIll 113 ■ 

This resemblance motivates us to study the gauged 
NJL model with extra dimensions. In this paper, we 
therefore focus on the dynamical chiral symmetry break- 
ing (DxSB) in the gauged NJL model with extra dimen- 
sions. The Lagrangian of the gauged NJL model with 
extra dimensions is given by 

G 



G r - 2 - ■ 



-itr(FM^F^^^) 



(II., 



where M,N = 0, 1,2,3, 5,- • • ,D and F^ is the chiral- 
ity matrix in D dimensions, r* {i = 1,2,3) the Pauli 



matrices, and G the four-fermion coupling. The gauge 
group is SU{Nc) and Fmn denotes the field strength. 
The gauge coupling is assumed to be on the UVFP. For 
simplicity, we take the number of flavor as Nj — 2, i.e., 
"0 — (V'lti V'd)"^ in tbc flavor space. The number of dimen- 
sions D is assumed to be even, D = 6, 8, 10, • • • , so as to 
introduce chiral fermions in the bulk. Extra S{= D — A) 
spatial dimensions are compactified at a TeV-scale R~^. 
It is convenient to rewrite the Lagrangian Eq. (jll.8|) by 
using auxiliary fields tr, tt^, 



C = ^ip^-^P{a + irAT'7T,)^-^{a^+Trf 



-^moa-itr(FM^F*^^) 



(II.9) 



In order to see the equivalence of Eq. (jII.Sp and Eq. (jII.Qp , 
we just need to eliminate the auxiliary fields a and tt^ 
through their Euler equations 



G 



■ipip + mo, t:. 



G 



ipiTAT'ip, (11.10) 



in Eq. (|II.9|) . The vacuum expectation value (VEV) of a 



G 

(o-) =-^(V'0) +"^o 



(11.11) 



is proportional to the chiral condensate in the chiral limit 
mo = 0. Note that the model with mo = possesses 
global SU{2)+ X S'C/(2)_ chiral symmetry. 

The bare mass term ttiq will be taken to be zero in the 
analysis of the D^SB in this paper. 



III. PHASE STRUCTURE 

A. SD equation 

We express the bulk fcrmion propagator as 

^S-\p)=A{-p')[^-J:{-p% (III.l) 

Nonvanishing S 7^ implies D^SB, which can be inves- 
tigated through the (improved) ladder Schwinger-Dyson 
(SD) equation. We need to take into account the run- 
ning of the bulk gauge coupling in the SD equation. In 
the SD equation there are three different momenta, (Eu- 
clidean) square of which are x = —p^, y = —9^ and 
z = —{p — q)^ for the external and the loop momenta 
of the fermion and the gauge boson momentum, respec- 
tively. Hence there exist various ways to incorporate 
the running effects. The simplest one [43 is to take 
5^4+5)13 (m) -^ ff^4+5)£)(™ax(^' y))' which has been widely 
used for four-dimensional models in Landau gauge where 



the SD equation dictates A{—p'^) = 1, so that it is consis- 
tent with the vector Ward-Takahashi (WT) identity. In 
the case of extra dimensions we also used this in Ref . [13 ■ 
However, this is not consistent with the chiral WT iden- 
tity |46|. and hence we adopt as in Ref. 11] an ansatz 
of Ref. |46l | , namely the gauge boson momentum is iden- 
tified as the renormalization scale of the gauge coupling 
strength. 



9{4+S)Dif^) ^ .9(4+<5)d(^) — 



g (m = V^) 

2-5/2 ■■ 



(111.2) 



where g is the dimensionless bulk gauge coupling. We 
assume here that the dimensionless bulk gauge coupling 
is on its UVFP, g = g*. The assumption is justified in 
case that the cutoff A is very large and the UVFP really 
exists. The SD equation is then given by 



A{x) 



A^ 



y 



D/2-1 



A{y) 



min(a;. y) 



X .L '^^A2(y)y + B2(y)[i^ax(x,y)]^/2-i 







nD-l}iD-4) 



(IH.3) 



B{x) = cr+{D-l+^)KD 
1 







^A^{y)y + B^{y) 



[max(x,y)]-°/^ 



(IH.4) 



where we have introduced ultraviolet cutoff A and the 
gauge fixing parameter is denoted as ^. B{x) is defined 
by 



B{x) =A{x)Y.{x). 
We denote the gauge coupling strength by kd. 



(IH.S) 



(in.e) 



with Cf being the quadratic Casimir of the fundamental 
representation. 



Cf^ 



N' 



1 



2Nr. 



(in.7) 



It is understood that a stands for the VEV of the aux- 
iliary field (ct) in Eq. Ijlll.41) . Hereafter, we will adopt 
this shorthand notation without discriminating the aux- 
iliary fields (a, tt) from their VEVs {{a), (tt)). When 
G — and hence a — mp, the SD equations Eqs. (|III.3|) 
and (|III.4|) are reduced to those of the bulk gauge theory 
without four-fermion interactions in a form given in (ll| . 
We choose the gauge fixing parameter ^ as |ll| 

{D-l){D-i) 



D 



(IIU 



With this choice of gauge fixing parameter, the fermion 
wave function factor becomes trivial, i.e., A{x) = 1 and 



B{x) = S(x). The SD equation then reads 

E(x) 



<J H ;:^ KD / dy- 



D 



y + E2(y) 



1 



The VEV of a is given by Eq. (jILlip . We thus obtain 



cr = ?7lo 



AD-2 



dxx^/2 



-1 S(x) 



-, (III.IO) 



/o a; + S2(a;)' 

where we defined the dimensionless four-fermion couphng 
constant g as 

and used 



[Nf = 2), (III.ll) 



{^^) = -NfTrSip) 






(III.12) 



The SD equation (|III.9p for the mass function is equiv- 
alent to the differential equation 

d^ ^, . D d 



^ X:2 ^(2^) + ^2;— I](a;) 



dx"^ 



2 dx 

2{D-l){D-2) x^{x) 



-kd- 



D '"^x + Y.^ 

with the infrared boundary condition (IRBC) 



0,(111.13) 



^C/2 "S(2;) 

dx 



0, 



(III.14) 



x=0 



and the ultraviolet boundary condition (UVBC) 
d ' 



{D/2-l)+x 



dx 



nx] 



x=A^ 



The derivative of Eq. Ijlll.91) at the cutoff scale A^ is given 
by 



x—E(x) 
dx 



2{D - 1)(D - 2) 



x=K^ 



DKD-2 ->^D 

X 1^ d.TX^/2-1 ^W 

/o x + S2(a;) 



From Eq. piI.10|l and Eq. piI.16|l we see 

D 9 d 

a = rrip, ; — x — i^yx) 

' 2{D - 1){D ~ 2) KD dx ^ ' 

The UVBC then reads 

D ( 



(Z?/2-l)+ 1 + 



A{D-1)kdJ dx 



= A2 



nx) 



(111.16) 



(III.17) 



=A2 



= (L'/2-l)mo. (III.18) 

For kd — 0, the SD equation reduces to the gap equa- 
tion of the pure N JL model. See the appendix of Ref. [13 
for the analysis of the D{> 4) dimensional pure NJL 
model. 



B. Bifurcation technique 

Even if we start with the chiral limit ttiq = 0, we ex- 
pect the dynamical chiral symmetry breaking S ^ takes 
place if the NJL coupling exceeds certain critical value. 
In order to fully analyze the behavior of the dynamical 
chiral phase transition in this model, we need to solve 
the nonlinear SD equation Eq. pil.9|l as it stands. Such 
a task turns out to be almost impossible within analyt- 
ical methods, however. Instead of solving the nonlinear 
SD e quat ion, here we first adopt the bifurcation tech- 
nique |4J| , restricting the solution much smaller than the 
cutoff scale S ^ A. The bifurcation technique is justi- 
fied at least for the determination of the critical coupling 
and the leading asymptotic behavior of the mass func- 
tion 'E{x). In this section, we employ the bifurcation 
technique for the sake of simplicity to avoid unnecessar- 
ily complicated expressions in discussing the decay con- 
stant, effective Yukawa coupling, and the propagator of 
the composite scalar, which are determined only through 
the dominant term of T.(x). On the other hand, the sub- 
dominant term is relevant to the scaling relation and the 
lowest interaction term of the effective potential. The bi- 
furcation method, however, does not always lead to the 
correct sub-dominant term. Thus we will perform more 
sophisticated analysis later in Section IIVI to determine 
the sub-dominant term of Y,{x). 

In the bifurcation technique, the S^ terms in the de- 
nominators of Eq. (|III.9|I and Eq. IJIII.10|) are replaced by 
an infrared cutoff Eq 

l^[x) = a -\ — kd I dy 



D 



^2 [max(a;,y)]^/2 1' 

(III.19) 



{D/2 - l)cr. (III.15) with 



JyD-2 



.A= 

/ dxx"/^-^J:{x), mo = 0. (III.20) 



We also assume So ^ A and Sq ~ S(Eo)- 

The integral equation Eq. I|III.19(I is equivalent with a 
set of the linear differential equation 

^2(£-2l(£^,,^,)^ 0,(111.21) 



D 



and boundary conditions 

^j:ix) 

ax 



= 0, 



X — Zjn 



and 



D Q \ d 
i:'/2 - 1 + 1 + -rrr^ -— I X- 



4:{D-1)kdJ dx 



S(x) 



(III.22) 



= 0. 

x=A^ 

(III.23) 



The solution of Eq. Ijlll.21|) behaves differently for 
i^D < '^S'* and for kd > k'j^^^. Here critical value k|5'* is 
obtained as pjj 



D D~2 
32D- 1' 



(III.24) 



For the sub-critical n^ < k™*i the solution of Eq. H111.21f) 
is given by a power-damping form, 



E(a;) =ciEo \ ^2 



-§(l-c^) 



dil^Q 1 2 



-^(l+UJ) 



X 

(III.25) 
where ci and di are real constants. We also define v and 

=3 



uj as 



and 



.4-. 



u;=J1-kd/k<=^'\ 



(III.26) 



(III.27) 



respectively. For the super-critical kd > K/?''; on the 
other hand, the solution starts to oscillate, 



T,{x) = ciSo J 2 
^1 



-^{l-iui) 



h.C, 



I X \ ^ I 1 X 

= 2|ci|So(^j sin ( -zyJ)ln^+. , , 

(III.28) 

with ci being a complex constant. The parameter Co and 
the angle 9 are defined as 



/ T 

y kd/k 



.crit 
D 



1, 



and 



2ie 



-ci/c 



1) 



(III.29) 



(III.30) 



respectively. 

Let us first discuss the case with kd < k™'. The IRBC 
Eq. 1III.22P leads to 



di 



1-LU 



Ci. 



(III.31) 



If we normalize the mass function such as S(I]q) — Sq, 
we can determine concretely ci and di. It is convenient 
to rewrite the UVBC Eq. (|III.23p in the following form. 



1 — uj"^ V I dx 



E(x) 



0. (III.32) 



x=A^ 



We use notations similar to those in Ref. l37 



Combining Eq. HIII.25|) and Eq. (|III.32ll , we obtain 



S|\ 



(1- w)ci -ig ~ iy{l + oj^ 



^AV {l+uj)di4:g-i^{l-uj)^' 

Eq. I|III.31(I then leads to the scaling relation^ 



^0 
A2 



4g- v{l +Lu)^ 



(III.33) 



(III.34) 



The nontrivial solution So 7^ exists only when the NJL 
coupling exceeds the critical line. 



r'* = ^(l+c.)^ 



(III. 35) 



See Figure ^ for the phase diagram in the n^-g plane. 
Note that the negative sign in Eq. (|III.31|I is essentially 
important. If the sign of di were identical to ci, we would 
not obtain positive value in Eq. Ijlll.33|) for g > g"^"^, 
indicating the instability of the vacuum. We will be back 
to this problem later in Sec. 

From Eq. Ijlll.25|) and Eq. HIII.31|) . we also obtain the 
asymptotic form of the nontrivial solution, 



■^1 



\{l-uj) 



for X > S^. (III.36) 



culated from Eq. (|III.15(I . We obtain 



We next consider the condensate a, which can be cal- 



(III.37) 



1+^ ^ /A 



2\ -■Ki-'^) 



neglecting the sub-dominant term. The asymptotic solu- 
tion Eq. Ijlll.36|) can be written in terms of a, 



2 /.x^-f(l-) 

^ ^ l+LO \A2. 



(III.38) 



The mass function S is thus proportional to the order 
parameter a. We will use Eq. HIII.38|) later in Sec. IIII E)l 
to calculate the Yukawa vertex of the composite scalar 
field a. 

We next turn to the case with kd > k^'*. From the 
IRBC Eq. (|III.22p we find 



2id l + ii^ 
e = —. 

1 — ILO 

On the other hand, the UVBC Eq. (|III.23|I gives 

/I A^ 

sin -z/o) ln—^ + + e'] =0, 

\2 Eq 



(III. 39) 



(III.40) 



* The scaUng relation Eq. iIII.341 is not correct near the pure NJL 
limit. Actually, it fails to reproduce t he s caling behavior of the 
D(> 4) dimensional pure NJL model ll(l in the Kjy — > limit. 
See Section llVl for more detailed analysis of the scaling behavior. 
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Broken Phase 




K 



FIG. 1: The critical line in {KD,g) plane. The critical line 
is given by ff'^"' ^ u{l +ujf/4 with uj = ^1 - kd /«:§'' and 
ly = D/2 — 1 for < kd < k^'*- The chiral symmetry is 
unbroken below the critical line. For k_d > k^'* the broken 
phase is realized irrespective of the size of g. 



Recall that the D^SB always takes place regardless 
the NJL coupling g for kjj > k™*. There does not exist 
critical NJL coupling with this kd- What does happen 
in the continuum limit for ko > k^'* then? Instead of 
g, one may think about tuning kd to its critical value 
K™*. This procedure cannot be done, however, since kd 
is essentially given by the UVFP of the gauge coupling 
strength in Eq. 



(jIII.6|) . The value of kd should be fixed 
once the discrete model parameters such as Nc and Nf 
are determined. We thus conclude that there is no fi- 
nite continuum limit for ko > k|5'' at least within the 
ladder approximation. We probably need to introduce 
other higher dimensional interactions in order to define 
the continuum limit with ko > k^'*. 

We also note that kb = k^'* cannot be achieved for 
D = 6,8 within the approximations we used.^ We thus 
concentrate upon the kd < nf,^^ case hereafter. 



with 6' being given by 



.-,^2 



^2i(s+e') ^ 4ff - v{l + iu) 
Ag — 1^(1 — iCb) 



(III.41) 



We note that nontrivial solution So 7^ exists in 
Eq. (|III.40p . irrespective of the size of the NJL coupling. 
This could be easily understood in an analytical manner 
by looking at the region a) ^ 1, where we can use an 
approximation 



1 

-r( 

VUJ 



V - Ag 



(III.42) 



So 


= exp 
= exp 


^ TT-{6 + 0') 


A 


^2 TT N 




yi^-Ag vCb) 



Eq. (|III.40|I then leads to a scaling relation (zero-node 
solution) 



(III.43) 



For g — vjA we find Sq/A — exp(a)/2) exp[— 7r/(2z^w)]. 

We next consider the continuum limit (A — > 00) in 
the present model. As we stated before, we simply ad- 
mit the existence of the UVFP in the gauge coupling 
strength in the present analysis. One might thus think 
that we should be able to obtain a finite theory even in 
the A — > 00 limit thanks to the UVFP. This is a non- 
trivial problem, however, in the D^SB vacuum. As we 
see in Eq. (jIII.34|) and Eq. (|III.43p , the dynamical mass of 
the fermion (So) is proportional to the cutoff A and thus 
diverges in the A — > 00 limit if both g and ku are fixed. 
In other words, in order to keep So finite, the NJL cou- 
pling g needs to approach its critical value Eq. (|III.35p in 
the A ^ 00 limit. This procedure, known as the Wilso- 
nian renormalization, allows us to define the continuum 
limit of the gauged NJL model. Hereafter we use the 
word "renormalization" in this sense. 



C. Decay constant 

There appear massless Nambu-Goldstone (NG) fields 
in the D^SB vacuum. The interaction of the NG field is 
then described by the low energy theorem and its decay 
constant. In this subsection, we estimate the size of the 
decay constant by using an approximation proposed by 
Pagels and Stokar (PS) J43|. The decay constant can be 
written solely in terms of the mass function S within the 
PS approximation. We know that the PS approximation 
works reasonably well in the D^SB of four-dimensional 
QCD. 

The D-dimcnsional generalization of the PS formula is 
given by 



Fl = 2^1'^Vl 



NDA 



cixX^/2 



-1S^(^)-§SS^(^) 



(x-fS2)2 

(III.44) 
where F-^ denotes the decay constant in D dimensions 



\Jm{^)VM)) = -igM^7r4 



(III.45) 



The four-dimensional decay constant /^ is obtained 
through the matching condition 



J TT 



{2t:RY 
2" 



Ft 



(III.46) 



It is straightforward to estimate the decay constant 
by using the asymptotic behavior of the mass function 
Eq. 1II06II . 



F^^S 



2+2i/(l-c^) 




dxx 



uu) — 2 



(III.47) 



For _D = 10 it happens kjj = k'S'*, for instance, for Nc = 81 and 



TV 



/ 



20. 



Note that F-^ diverges in the A — > oo hmit for 



1 

UJ> -. 
V 



(III.48) 



The condition Eq. (|III.48p can be rewritten in terms of 



>^D < ^T, 



(III.49) 



with the marginal triviahty (MT) point k^"^ being de- 
fined as 



'^D^ - ( 1 lyi ] '^D 



(III.50) 



The NG field interaction is suppressed by 1/F.^ accord- 
ing to the low energy theorem. For kd < k^'^ the NG 
field is therefore decoupled from the rest of system and 
becomes trivial in the continuum limit (A -^ oo). 

On the other hand, surprisingly enough, the decay con- 
stant Fjr remains finite for 



MT 
ft, 25 



< KD < K™ , 



(III.51) 



even in the continuum limit A — > cx). Once the condition 
Eq. pil.51|) is satisfied, the NG field enjoys non-trivial 
interactions in the "renormalized" theory of the higher 
dimensional gauged NJL model. 

The MT point k^'^ is indicated by a cross in the phase 
diagram Figure ^ At this point the decay constant F^^ 
diverges as F^ ^ In A, while F^ ^ A^^'"^"^) for < k^ < 



where Tg {q + p, q) is the Yukawa type vertex of the 
auxiliary field, and q and p are the "0 and a mo- 
menta, respectively. The basic difficulty in the calcu- 
lation of Eq. (|III.52|I arises from our lack of knowledge 
of Ts{q + p,q). In Appendix A, we calculate analytically 
Ts{q + p,q) with certain approximations. As for the ef- 
fective Yukawa vertex at p = 0, on the other hand, we 
easily find 



da 



1+LO 



-q 

A2 



2\ -§(1-'^) 



(III.53) 
where we used Eq. (|III.38p . For the relation between E 
and Fs, see, e.g., Ref. 50]. 

ATW showed that Eq. 111152^ can be approximated 
by 



zD-\p) = 



d^'q 

i(27r)^ 



tr 



r,(-q2)ir,(™g2)-i 



-constant. 



(in.54) 



within reasonable assumptions. Note that Eq. Ijlll.54|l 
contains Yukawa vertex only at p = 0, 



F,(-g2)=r,(g,(j), 



(in.55) 



and we are thus able to use Eq. (JIII.53P in the evaluation 
of the a propagator. 

After the Wick rotation and performing the angular 
integrals, we obtain 



In D —f 4 we find kj^^ -^ 0, namely the MT point iD^^(n) — iD^^(O) 



coincides with the pure NJL point, where it is well-known 
F^ diverges logarithmically. On the other hand, in the 
D ^ oo limit k^^ -^ ^^S'*' ^^"^ hence the region of the 
finite -FV is squeezed out. 



D. Yukawa coupling and scalar propagator 



-2°/^NcNf 



d^q 

i{2tt)D 



rii-q') 



2^P-<1-P 



q^ip-qf 



where we defined 



In addition to the NG bosons, a composite scalar boson 
also appears in the spectrum of the NJL type models. As 
we will show later, the auxiliary scalar field a of Eq. I|II.9|) 
actually acquires non-trivial propagator at the loop level. 
In this subsection, we adopt a simple method proposed by 
Appelquist, Terning, and Wijewardhana (ATW) 48] to 
calculate roughly the properties of the composite scalar. 
We also describe more sophisticated method based on 
Ref. [43 in Appendix A. These results qualitatively agree 
with each other. 

The auxiliary field a propagator at \p'^\ ^ Sg is given 
by 






i{2n)^ 



tr 



Fs(g+p,g)-7 



-^ 



-t-constant. 



(III.52) 



-P 



The kernel K is given by 



y 



-q 



K{x,y) = 
2 min(a:, y) 
D max(x, y) 

X 



F[2,2-D/2,D/2 + l; 



^{x,y) 



max(a;, y) 



max(x, y) 
For even dimensions we obtain 

.5/2+1 



F ( 1, 2 - D/2, D/2; ™''^,^'^^, ) . (III.57) 
' max(x,y) ' 



K{x,y) = -Y. 



e=o 



(-S/2-l), ^yY 
{S/2 + 2)e \x 



(!) "<--") 



fw^ei»<»-^''<"-' 
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where {a)g is defined as 

{a)i ^{a + i- l)(a + £ - 2) ■ ■ ■ {a + l)a (III.59) 



and (a)o = 1. Substituting Eq. (|III.53p for Eq. (|III.56p . 
we obtain tlie composite scalar propagator for uj ^ 1/v 
{no ^ ^T)-- 



tD-\p) - iD-\0) = 2''/^N,Nfn 



NDA 



l+OJ 



xA 



D-2 



c.j{,r*"j:cj-y 



e=i 



vav 



where 



a 



and 



S/2+1 



1/1,: ^ — ' 



(-(5/2 -1)£ 2vuj 
(^2 + 2)^ (i/a;)2-£2 



,(111.60) 



(III.61) 



G 



_(-V2-l), 1 
(J/2 + 2)£ vuj-V 



> 1). (III.62) 



At the marginal triviality point oj — 1/v {kq — i^^^): 
we find 



iD-\p) - tD-\0) = 2'"''N,Nfn 



NDA 



l+LO 



xA 



D-2 



Q 



"(^) 



In 



A^ 

X 

S/2+1 



<^1 V \2 



^^■ii) 



where 



and 



Co^-- 



e.=2 



D-2 



D 



(III.63) 



(III.64) 



r' - -^ - ± - 'v' (-V2-l)£ 2vuj 
'"2 i^c^ ^ (V2 + 2)£ (:.a;)2-£2- 

(III.65) 
We next consider the A ^ oo limit. It is convenient to 
renormalize the auxiliary field <t 



OR = Z-i/V, 



(III. 



in such a limit. We define the renormalization constant 
Zc so as to keep the renormalized propagator finite, 



z- 



2i/(l-w) 



KY'"-'' /A2 
ln( — 



2(1.-1) 

A\2(''-i) 



(III.67) 



The renormalized propagator is then given by 



2D/^NcNfn 



NDA 



,2iy 



\^2 



^OM^'' ( ^ ) , 



,2iy 






(KD = K^ ), 

(0 < /.z? < ^D, 



(III.68) 



in the A — > oo limit. 

What does happen for the Yukawa type vertex, then? 
From Eq. HIII.53() and Eq. (|111.67|l we obtain 

Vf\~q^)^Z\I^V.[-^^)^ 



(III.69) 

in the A ^ cx) limit. We note that the composite 
scalar thus decouples from the rest of the system for 
>^D < '*!) in the continuum limit. On the other hand, 
for K^""" < kd < Kg'*, the composite scalar interacts with 
fermions in a non-trivial manner through the Yukawa ver- 
tex Eq. (|III.69p . These behaviors are consistent with our 
analysis of the NG field (decay constant) in the previous 
subsection. 



IV. MORE ABOUT THE MASS FUNCTION 

For the determination of the precise scaling behavior 
around the critical line, it is not enough to know the dom- 
inant asymptotic behavior of mass function Eq. (|III.36p . 
Actually, lacking the sub-dominant term, the bifurcation 
method does not reproduce exactly the same scaling re- 
lation (gap equation) as that of the pure NJL model in 
the vanishing gauge coupling limit even in four dimen- 
sions |51|. We also need information of sub-dominant 
terms. The sub-dominant terms are also required in the 
calculation of the effective potential. Unfortunately, the 
bifurcation technique is not enough for such a purpose. 
In this section, we evaluate the sub-dominant terms first 
in the well-known linearized SD equation, and next in 
the power expansion method. 

It is known that the linearized approximation works 
well in four dimensions. By using the linearized ansatz, 
we can solve analytically the ladder SD equation. We ob- 
tain the coefficient of the sub-dominant term as well as 
its exponent. The power expansion method reflects the 
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correct asymptotic behavior of the nonhnear SD equa- 
tion. However, we can determine only the exponent of 
the sub-dominant term in the method. The coefficient 
of the sub-dominant term is left as an unknown param- 
eter. In a certain region of the gauge coupling, it turns 
out that the exponent of the sub-dominant term is dif- 
ferent depending on the methods. We thus describe both 
results. 



A. Linearized approximation 

We first consider the linearized SD equation [32. |43 | , 
in which the Y? term in the denominator is replaced by 
a constant Sq = Y?{x — 0). This approximation works 
well in the determination of the dominant asymptotic so- 
lution at a; ^ Sq. It is also known that the scaling behav- 
ior determined from the linearized SD equation closely 
approximates to the numerical results of the non-linear 
SD equation in the case of the four-dimensional gauged 
NJL model. 

The linearized SD equation is given by 

v( ^ ^ 4(J-1) f'^' y^l^-'nv) 



D 



y + s 



c{x,y)]Dn 



(IV.l) 



(T = Too 



and a subsidiary condition. 



dxx 



^/^-'^%,(IV.2) 



So = I](a; = 0). 



(IV.3) 



The integral equation Eq. (|IV.1|I is equivalently rewrit- 
ten in terms of a set of differential equation and boundary 
conditions, 



dx"^ 



x'^-—i:+{v+i)x—i: 

ax 



1/2(1 -u;2) xY. 



-^^Yix) 
ax 



x + Y'' 



= 0, 



d 

V + X—- 

ax 



S(x) 



x=a 



= A2 



va, 



0, (IV.4) 

(IV.5) 
(IV.6) 



where v and w are defined in Eq. HIII.26(I and Eq. (|111.27|l . 
The UVBC Eq. IJIV.6|) can also be expressed as 

4 



1 



i\^d_ 

\ — up- V j dx 



Y{x) 



vma. 



x=A^ 



(IV.7) 

Note that the differential equation (|IV.4|I possesses 
three regular singular points at a: = 0, — I]q,oo. The 
solution is then expressed by using the Gauss hypergeo- 
metric functions. It is now easy to find the solution of 



Eq. (|IV.4|) satisfying the subsidiary condition Eq. IJIV.31) 
and the IRBC Eq. PV3|| . 



E(x) = So^ ( ^(1 + c^), ^(1 - ^), ^ + 1; -^ 



where F is the Gauss hypergeometric function 



F(a,/3,7;z) 



v^ {a)n{l3)n z 

n=0 



(7). 



(IV.8) 



(IV.9) 



with (a)„, (/3)„, (7)„ being defined in Eq. (|III.59|I . We 
assumed here K£, < k^'* and hence a real and positive lu. 
The UVBC Eq. (|r\r7|l determines Sq and thus the scal- 
ing relation. We need to know the behavior of the mass 
function Eq. IJIV.8|) in the ultraviolet region a; ^ Sg. 
For such a purpose, it is convenient to use a well-known 
formula |52 |. 



F(a,/3,7;z) = 
r(7)r(/3-a) 
r(/3)r(7-a) 

r(7)r(a - (3) 
r(a)r(7 - /3) 



{-z)-"F{a,a--f + l,a-(3- 



z 

z 
(IV.IO) 

Combining Eq. HIV.9|I and Eq. (|lV.10f) . we expand the 
solution Eq. HIV.8|) for a; > Sg, 



S(x) = Eo ( :^ 22^'" 

^0/ „=i 



+So 



X 
^0 



Ki+'^) 



J2d„ 



X 



-n+l 



(IV. 11) 



with c„ and (i„ being given by 

(-1)"-! Tiiy + l)T{iyLu) 



(n-l)!r(|(l + c.))r(f(l+c.) + l) 
(§(1 --))„-! H(l +-))„-! 



d„ = 



i-iy 



T{iy+l)T{-i'Uj) 



(IV. 12) 



(n-l)!r(|(l-c.))r(f(l-c.) + l) 



(f(l+^))n-l(-f(l~^)). 

(^'^ + 1)„-1 



(IV. 13) 



The ci term in Eq. IJIV.11|) gives the leading contribu- 
tion to the mass function T,{x) in the asymptotic region 
a: ^ Eq for kd < k™'- The next-to- leading contribu- 
tion, on the other hand, comes from C2 or di , depending 
on the value of kd . For k^'^ < ke, < k^'* di term gives 
the next-to-leading contribution, while C2 term becomes 

"^ the 



..MT 



the next-to-leading for kd < k^'^ ■ At kd 

power of C2 and di terms goes same in the expression of 

Eq. HlV.llfl . Although the coefficients C2 oc l/(— i/w -I- 1) 
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and di ex r{—Piv) diverge at kd = '^15'^, i-e-, w — 1/j^, 
they cancel out each other, so that the logarithmic term 
appears, 



E{x) 



ClSo ( T^ 



X 
2-1 



-|(''-1) 



I/' - i / X 



-^(''+1) 



In 



y2 
^0 



{KD = n^j^'). (IV. 14) 



In order to determine the scaling relation near the crit- 
ical line, we compare the leading (ci) and the next-to- 
leading (di or ci depending on kd) terms in the UVBC 
Eq. pV!7|) . For mo = 0, wc find 



So 

A 



4w di 



,crit 



1-^— ), (IV.15) 
9 



for K^^ < Kd < n'g , and 



So 
A 



iy{l-Lu'^)ci 



C2 



1- 



(IV. 16) 



for Kd < Kd'^ , where we used g ~ 5'^'^'*. Note here that 
the signs of ci and di (02) are opposite to each other for 
K^'^ < KD (kd < K^^)- We thus find 



So 

A 



, (k^ < kd < Kd ), 



iO<KD< K^-^) 



(IV.17) 
Note here that Eq. IJIV.17P in kd -^ limit agrees with 
the scaling relation obtained in the £>(> 4) dimensional 
pure NJL model lOj. 

The chiral condensate is calculated from Eq. (jIILlSp . 






U(l-Lj) 



l^il+i^) 



(IV.18) 



for K^T < K£, < Kg", 



and 



(1+^) ^ /So 



u{l-uj) 



l+UJ 



I) ^^^° (t 



i^(l-w)+2 



(IV. 19) 



for < KD < K^^- 

Further details of our analysis on the linearized SD 
equation are given in Appendix IbI 



B. Power expansion method 

We next try to investigate the sub-dominant asymp- 
totic solution directly from the non-linear SD equation 
Eq. Ijlll.9|) without using the linearizing ansatz of the pre- 
vious subsection. 43] The non-linear differential equation 
Eq. HIII. 13(1 and the boundary conditions Eqs. HIII. 14(1 
and piI.15(l can be rewritten as 



.d^E 



^S:3 + (^ + i)^:7r 



dS 1/2(1 -cj2) xT, 



dx^ 



X' ''— S(a;) 
dx 



dx 
= 0, 



a;-^S2 



0,(IV.20) 
(IV.21) 



x=0 

4 g 



1 — UJ-^ V 



dx 



S(x) 



vmo. (IV.22) 



x=A'^ 



For X ^ Y? {x) the last term in the non- linear differential 
equation Eq. (|IV.20|) can be expanded as 



xY, 



a; + S2 



S3 £5 
= S + — 



(IV.23) 



Eq. (|IV.23p allows us formally to expand the solution E 
of Eq. JiV.2()ll . 



S(x) =Si(x)-f E2(a;) + ' 



(IV. 24) 



where Si (a;) is the solution of a homogeneous linear dif- 
ferential equation 



■ 2 d2 ^ 1.2(1 _ ^2) 



Si(a;) =0, 



(IV.25) 
while E2 (a;) is the inhomogeneous part of the solution of 



dx dx 4 

_ uHi~Lj^)j:l(x) 



Y2{X) 

(IV.26) 



The higher terms in the expansion Eq. (|IV.24I) are given 
in a similar manner to Eq. pV.26(l just by replacing the 

R.H.S. to S?E2, SiEi, E^, Ef, -^ 

What does the expansion Eq. (|IV.24|1 mean then? The 
solution of Eq. I|IV.25|) is given by 



Ei(x) =: ciEo I 2 



-^i^ii-u,) 



diEo 



E2 



-^iy(l+Lo) 



(IV.27) 

with Eq being a constant with a mass dimension. We 
thus obtain 



S2(X) = 
C2S0 



E2 



-fi.(l-cj)-l 



CoSo 



E2 



-rf2So ( ^ 



|i.(l+w)-l 



\u+^i^uj~l 



lu—irUU—l 



d^So(^) , (IV.28) 
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with C2 oc cf, C2 ex c^di, ^2 oc df, ^2 ex cid^. Note that 
the exponent of the power damping of 22(2;) is steeper 
than that of T,i{x) for x 3> Sq. It is also obvious that the 
higher terms in the expansion Eq. IJIV.24|) falls quicker 
even than S2(a;) for a; :» Eq. The expansion Eq. HIV.24|) 
can therefore be understood as an "power expansion" in 
this sense. 

Since we are interested only in the dominant term (ci) 
and in the sub-dominant term {di or C2 depending on 
kd), we restrict our analysis to these three terms. The 
power damping behavior of ci and di terms agree with 
those in the linearized SD equation. We note, however, 
the exponent of the power damping of the C2 term differs 
from the the result of the linearized SD equation. Actu- 
ally, the C2 term becomes sub-dominant for kjj < k^, 



.PE _ 



-K-^ 



We also note 



.PE 



< K 



MT 
D J 



(IV.29) 



(IV.30) 



for D> A. 

We first consider the case kd > k^, where the di term 
gives the sub-dominant contribution. Since the power ex- 
pansion cannot be adopted in the infrared region, we are 
not able to use the IRBC Eq. HIV.2ip to determine the 
coefficient ci and di. If we assume that di possesses 
opposite sign to ci, we find a scaling relation similar to 
Eq. IjIV.lSp. F or di/a > 0, on the other hand, the R.H.S. 
of Eq. (|IV.15p becomes negative. Eq. IJIV.15(I thus can- 
not be understood as a scaling relation. This problem 
may be related with the vacuum instability which will be 
discussed in Sec0 Hereafter we assume the result in the 
linearized SD equation, di/ci < 0, remains to be valid 
even in the power expansion method for kd > k^ . 

We next turn to the kd < k^ case, where the C2 term 
gives the sub-dominant contribution. From Eq. ()1V.26[I . 
we find 



C2 



i.2(l-c^2) 



4(1/ + 1 - iyLL!){2iyLU - J/ - 1) 



-ll 



(IV.31) 



and thus C2/C1 < 0. The scaling relation is then given by 



As in the linearized approximation, the coefficient C2 in 
Eq. HIV.31|) diverges at kjj = k^, i-e, w = 1/2 -I- 1/(2:/). 
This implies that the logarithmic term appears at the 
point. We can easily confirm that S2(a^) given by 



^2{X) 



(^-1)(3^ + 1) 3 



(^+1) 



ci^ol^j In^, 



(IV.34) 

satisfies Eq. (|IV.26|I . 

The chiral condensation is calculated from Eq. (|111.15|l . 



1 /A2 



A^ 
-{l-u;)diY.o ( ^ 



-f(i-t^) 



2\ -■1(1+'^) 



(IV.35) 



for K^ < Ku, while 



-(1 +uj)ciT.Q 



A2 



Ki-^) 



i{'-'-^l)''-''^{^. 



-(l-uj)-l 



(IV.36) 



for < KD < K^- 



V. EFFECTIVE POTENTIAL 

We here outline our method to calculate the effective 
potential V{<7) in the gauged NJL model. See Refs. pH 
i53, 54] for details. For simplicity we assume tt.^ = in this 
analysis. It is easy to restore the tt^ degree of freedom by 
using the chiral symmetry, i.e., a'^ —> a'^ + irf . 

We start with the partition function VF[J], 



W[J] = - In / [dV'^] [gauge] exp U / d^x{£ + J a) 



(V.l) 
with C being the Lagrangian of the auxiliary field descrip- 
tion of the gauged NJL model Eq. I|II.9|) . For a constant 
external field J, the partition function can be written as 



So 

A 



2v{l-Lj)+2 



v{i - ^2) ci ( ^_ .r 



4(1/ -|- 1 — vu) C2 \ g 

with this Kd , where we used g ~ g'^"* . We thus find 

1 

^crit \ 2u^ 



So 

A 



1- 



9 



9 



(IV.32) 
We thus f] 

/..PE ^ ,. ^ ,^crit\ 
{hjj < h,D < Kjj j. 



9 



crit \ 2i/(l~ij)+2 



1-2 ^ 0<Ki3<k£^. 

9 J 

(IV.33) 



We note 



W[J] 



dJ 



d xw{J). 



^^{J) = a, 



(V.2) 



(V.3) 



with a in the R.H.S. being understood as the VEV of cr. 
Solving J as a function of <t in Eq. (|V.3p . we define the 
effective potential V{<7), 



V{a) ^Ja-w{J), 



(V.4) 
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and thus obtain 



7 t^(7 

-j-V{a) = J, or V{a) = / dcrj. 



(V.5) 



We next describe how we calculate J in terms of a. 
Looking at Eq. (|II.9|I and Eq. I|V.1|I . we find the effect 
of constant J can be taken into account by replacing the 
bare mass mo, 



-^rriQ -^ -^mo + J. 
in the SD equation. From Eq. (|III.17|) we obtain 






-ma + a + 



1 — uj"^ v"^ dx 



x-—Yi{x) 



(V.6) 



(V.7) 



a;=A2 



The derivative at the cutoff ^I](a:;)| _.^ can be calcu- 
lated from the UVBC Eq. 1111.1511 . 



x-^E(x) 
dx 



x=K^ 



u [a - E(A2)] , (V.8) 



which leads to a compact expression for J, 



^^ ^ '[a-m^)]. (V.9) 



iV, 



■J = —mo + (7 + 



1 — uj-^ V 



We thus obtain the effective potential V{a) in terms of 
a and ^(A^): 



2^/^NrNfK^Vi 



cl^fl 



'NDA 



V{<j) 



1 m,o(T 1 CT^ 

+ ^;^i,/d.[a-E(A^)]. (V.IO) 
In the pure NJL limit kd = 0, i.e., oj = 1, we should 



return to Eq. pil.9(l . Plugging with T,{x) = a and 
Eqs. (|III.10p . (EH), we easily find 



'2^/2iV,iV^A^rj 



'NDA 



Via) 



1 mo<J 1 g'^ ^ I a 



1 moCT 1 / 1 1 \ ''' 



dz- 



z - 



1 a^ 



g A2 2 V.g !^/ A2 4(i/-l) A4 



+0(^1, («;z5 = 0). 



(v.ii) 



Since the sub-dominant term of ^(A^) with nonvanish- 
ing gauge coupling depends on the approximations, we 
calculate V{(t) separately. 



A. Linearized approximation 

In the calculation of the effective potential V{(t), we 
need to express ^(A^) in term of a. Let us start with 
the solution of the linearized SD equation for k^ > k^"^. 
Combining Eq. IjlV.lSp with 



E(A2) = ciSo 



we obtain 



A 



iy(i-w) 



diEo 



A 



1^(1+^) 



(V.12) 



E{A^) = ^—a- 



2to , , / 2c7^ a 
-diA I ^ 



■^(i+'^)+i 



1+Uj" ' 1 +U}""'' \1 +U) A 

(V.13) 
We are now able to calculate the effective potential from 



Eq. 1V21), Eq. SSlM and Eq. JSTm 



2^/27Ve7V/A^r!NDA ) Via) 



1 7710(7 1/1 1 \ O" 



5 A2 2 \g g''"^ J A^ 



Mj 



,(V.14) 



for K^^ < kd < K^ . Here the coefficient Ai is defined 
as 



Ai 



■di ' 



l-Cj2l+cJ I/(i^+l) """Vl+^/ 

(V.15) 
It is straightforward to perform similar analysis for < 
Kd < K^'^. We obtain 



(2^/2iV,A^/A^f7NDA) V(a) 
1 m-ocr 1 / 1 



5 A2 2\g g"^^J A2 

2+T77Tin 



--=(x 



(V.16) 



for < Kd < K^'^ . The coefficient A2 is defined as 



A. 






2(i^w - 1) i^(l - tj) + 2 ^Vl+^ 



(V.17) 

The calculation at kd — k^'^ is a little bit involved, 
since we need to take account of both di and C2. By using 
formulas given in Appendix lEl we obtain 

(2''^^ N,Nf A^'n^^A)'' Via) = 

1 mocr ,1/1 1 \ 0-2 

gClit J ^2 

,„^ 2+2 / ^ X 

•••,(V.18) 



g A2 2\g g^^y A^ 



-Mj) 



a 
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for K£, = K^ and the coefRcient A'2 is given by 



A' 



y + l 



2i/Ci 



i\ 2+1 



(V.19) 



Using Eqs. (|IV.12|I and (|IV.13|I . we find the coefficients 
Ai, A2 and A'2 are ah positive definite in the hnearized 
SD equation. The potential is therefore stabilized for 
large value of a even if we truncate the potential and 
neglect H terms in Eqs. (|V.14|I . (|V.16|) and (|V.18|) . 



B. Power expansion method 

The effective potential V{(j) can also be calculated in 
the power expansion method. We should keep it in mind, 
however, that the coefficients ci and di cannot be de- 
termined with the power expansion method, unlike the 
calculation in the linearized SD equation. 

We here summarize our results of the effective potential 
in the power expansion method: 



1 moCT 1/1 1 \ ^^ 

~5^V2~^ 2 \g^^) A2 



+ii(5) "^-"^^+---,(V.20) 



for K^ < kd < K™*i ^-iid 



h'^/^NcNfA^n 



NDA 



V{a) 



1 mocr 1/1 1 \ o" 



g A2 2 \g g"'^ J K^ 



-^ ^ + • • • , (V.21) 



for < K_D < K^^, and 



(2^/2^ciV/A^r!NDA) V{a) 



1 TTIqCT 1/1 1 \ <'' 



g A2 2\g g"'^J A^ 



+^2 ^ In (^ )+•••' (V-22) 



for kd — kW". The coefficients Ai, A2 and A!^ are defined 



by 
ii = 



■/(i+c^)+i 

■rfi 



l-Cj2l+cJ i^(i^+l) 



1 +tJ 



1 



A{2vijj-v-l) \\+ioJ ' 
4 



1 



41/ 



1/ + 1 V 3i/ + 1 



(V.23) 
(V.24) 

(V.25) 



respectively. 

We emphasize here, due to the lack of our knowledge 
of Ci and di , the coefficient Ai cannot be calculated in 
the power expansion method. On the other hand, we find 
A2 > for < kd < K^, independently of ci. 

We comment that the exponent of the lowest interac- 
tion term {Ai term) in Eq. (|V.14|I is the same as that 
of Ai in Eq. IJV.20|) . In passing, the bifurcation solution 
also leads to the same exponent of the lowest interaction 
term. Thus the lowest interaction term in the effective 
potential is insensitive to the approximations of the SD 
equation for k^'^ < kd < k^'*- ^°t^ that k^ < k^'^ . 
On the other hand, for < ku < k^"^ , the exponent of 
the lowest interaction term is different depending on the 
approximations. 



VI. NONTRIVIAL WINDOW 

In Sec. mil we found "renormalizability" or "nontrivi- 
ality" in the sense that the decay constant F^ and the 
renormalized Yukawa coupling Tl become finite in the 
continuum limit for k^"^ < kd < k'j^^^. 

Here we further study the renormalization and discuss 
the renormalization group (RG) flow, based on the effec- 
tive potential derived in Sec. The renormalization of 
the effective potential is performed in a way similar to the 
four-dimensional gauged NJL model [33 ■ This method 
is applicable not only in the broken phase, but also in 
the symmetric phase. Note that we already renormal- 
ized the a field in Eq. l|III.67p so that its kinetic term is 
kept finite. 

For K^"^ < Ku < K^'' we require the bare parameters 
(mo and g) in the effective potential to depend on the 
cutoff (mo (A) and g{A)) such that 



3(A) 



const. 



and 



1 



1 



5(A) g' 



crit 



A^""^ = const.. 



(VI.l) 



(VI.2) 



The effective potential is thus kept finite even in the con- 
tinuum limit A ^ CX3 (Wilsonian renormalization). We 
define the renormalized parameters m.^, 5^ as 






mo (A) 
.9(A) 



^iy{l+uj) 



(VI.3) 
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and 



1 



M 



2uu 



1 



A' 



2i>ci; 



9R{^^) 5"" / V 5(A) g"'* 

(VI.4) 
with /i being the renorniahzation scale. The effective 



potential Eq. IJV.14(I (linearizing approximation) is thus 
renormalized as 



(2^/2iV,iV^^^rjNDA) ' Vr{cJr) 



1 niRaR 1/1 1 \ <7f? 






+^l(^ 

M 



, (VI.5) 



for K^T < K^ < K^jit, 



and similarly is Eq. (|V.20|) (power 

Note that the "H " terms in 

are decoupled in the A 



expansion method) 

Eqs. HV.14|I and l)V.20ii are decoupiea m tne iv ^ oo 
limit. The a field has non-trivial self- interaction in the 
renormalized effective potential. Recall that Ai > 0. The 
potential is therefore stabilized and the VEV of aR re- 
mains finite for gR > g"^^. We emphasize that the form 
of the renormalized effective potential VR{aR) does not 
depend on the approximations, although the coefhcient 
of the interaction term does. 

From Eqs. (|V1.3|I and (|V1.4(I we obtain the beta func- 
tion of the dimensionless four-fcrmion coupling gR, 



Pign) = '^i"^9R 1 



gR 

^crit 



(VI.6) 



and the anomalous dimension of the fermion mass "fm , 






gR 



l~uj + 2uj 



gR 



g^ 



(VI. 7) 



for K 



MT 
D 



< Kd < k' 



crit 
D ■ 



They take the same form as 



those in the four-dimensional gauged NJL model with 
fixed gauge couphng [23 up to the factor i/ = D/2— 1 -^ 1 
(D — » 4). We note the beta function of the NJL coupling 
possesses a UVFP gR = 5"'* (= v{l + tj)2/4). At the 
fixed point of gR the anomalous dimension 7^ reads 



7m = iy{l + 10). 



(VL8) 



Let us turn to the region < kd < 1^^'^ ■ In this region 
the decay constant -FV diverges as we mentioned. It is 
thus expected that the same renorniahzation procedure 
as in the region k^"^ < kjj < k'j^^^ will break down. In 
fact, if we formally performed such a renorniahzation, we 



would get 



(2"/^NcNffi^flNDA) ' VRiaR) 



1 niRaR 



1 

gR 



gR M 

where we defined ttir and gR as 



mRi^i) ^+i^mo(A)^ 



f+i 



and 



1 



gRifj-) 



1 



gRifj-) g'' 



5(A) 



1 



(VI.9) 



(VI. 10) 



.9(A) g 



1 

crit 



3^)A^ (VLll) 



This renorniahzation would imply the beta function and 
the anomalous dimension: 



and 



P{9r) = 25fl 1 



7m(5fl) 



gR 

^crit 



(VI. 12) 



/3(ff. 



R) 



gR 
1^-1 + 2 



gR 



■crit 



g' 

The anomalous dimension 7,„ would read 

7m = i^ + 1 



(VI. 13) 



(VI. 14) 



at the UVFP of gR. 

However, the self-interaction term in VR^aR) 
Eq. (|VI.9p disappears at A — > cxd in our renormal- 
ization procedure both in the linearized approximation 
and in the power expansion method. Thus the renor- 
malized potential is not stabilized for gR > g"^"*^ and 
hence the renormalization breaks down. The Yukawa 
interaction does also vanish and hence interactions of a 
and TT are all trivial. 

Let us next consider the dynamical dimension of 
composite operators. Within the ladder approximation 



® At K£) = K,^"^ the renormalized effective potential would take 
the same form as Eq. iVI.9i . Although the definition of gn and 
rriji in Eqs. iVl.lOl and IVl.lli should be modified by loga- 
rithmic factors, the beta function and the anomalous dimension 
are unchanged from Eqs. IV1.12i and IV1.13i . If one took the 
limit Ko \ K^"^ of the expression in Eq. iVI.5l . one would find 
the pole in the coefficient Ai. This pole actually is logarithmic 
divergence In A when the sub-sub- leading term is properly incor- 
porated in the limit A —> 00 as in t he pure NJL limit of the 
four-dimensional gauged NJL model |37|. After the renormal- 
ization of the <T field Eq. IIII.67i . we find the interaction term 
vanishes in V[i{aji), consistently with Eq. iV1.9i . 
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adopted in this paper, various composite operators have 
the dynamical dimensions, 

dim{{i>ij)^) = 2dim(i/iV'), (VI.15) 

dim((ViV')'') == 4dim(VjV'), (VI.16) 

dim(5Af(ViV)9^(V'V')) = 2dim(V;V') + 2,(VI.17) 

with 

dim(V;V') = ^ - 1 - 7m. (VI.18) 

Using Eqs. ifVLS)) and (|VI.14|I . we find 




K 



dim((V'V')^) = 



D-2i^Lo, {k^'^ < KD < k'E''), 



D-2, 



(VI.19) 



dim((V'V')'') 



2D - AviO, (k^T < k^ < ^crit)^ 



2i?-4, (0<K^<A.MT), 



(VI.20) 



D + 2- 2vuj, 



(k^t < ^^ < ^git)^ 



(0 < ^^ < kMT) 



(VI.21) 



The four-ferniion operator {'ip'ijj)'^ is relevant in the whole 
region < kjj < kJj'*' which justifies the inclusion of 
the four-fermion interaction in the study of the phase 
structure of walking gauge theories. 

For the_ region Kj^'^ < kd < k™* both operators (V'V')^ 
and dM{4'i^)d^ i^i^) are irrelevant, i.e., dim((-0-0)^) > D 
and diin{{dM{i^^)d^ {jjjip)) > D. This is consistent with 
our renormalization procedure without introducing such 
operators. 



On the other hand, for < k/j < 



.MT 



the 



kinetic operator dM{''l''4')9'^ (ipip) becomes marginal, 
dim((9M(V''0)c^*'^('0V')) = -D, while the eight-fermion op- 
erator {'ipip)'^ is irrelevant, diin((?/j'0)'*) > D. This is an- 
other symptom of breakdown of our renormalization in 
this region. We may need to incorporate the marginal 
operator 97\f ('0V')9*^(V'?/') as a counter term, which is 
similar to the situation in the four-dimensional pure 
NJL model where both (9p(V'?A)9''(V''0) and {ipip)'^ are 
marginal l37l|. 

We comment on further symptom of breakdown of our 
renormalization procedure for < kd < i^^"^ ■ If we 
renormalized consistently the effective action in this re- 
gion, we would expect the renormalization of ipijj as 



mR{{'ipip)R) =r7io((V"0)o) 



FIG. 2: Phase diagram in {KD,gR) plane. The critical line 
g"^^ is the same as that in Fig. Each point on this line 
describes a different theory having different Nc,Nf. In the 
region (A) kd > Kd'* the dynamical mass diverges in the 
continuum limit, i.e., no finite theory exists. In the region 
(B) K^^ < Kd < tf)'* (nontrivial window), the dynamical 
mass can be made finite and the self-interaction of an and the 
Yukawa interaction remain nontrivial in the continuum limit. 
In the region (C) < kd < k^"^ the renormalized effective 
potential is not stabilized for the broken phase qr > g"^^ and 
hence the renormalization breaks down. In this region the 
self-interaction of (t_r and the Yukawa interaction vanish and 
are trivial. 



Note that mo ^ A"''- = A^'^^i, while Eq. (|III.12p yields 



{i4>-ip)o) oc 



/' 



dxx^ 



nx) 



^u(U 



(VI.23) 



where we used the asymptotic behavior of the dominant 
solution in Eq. HIII.36I) . Thus Eq. (|VI.22p diverges as 
mo{{tpilj)o) ^ A"'^"^. This can be understood by not- 
ing ('ijj'ip)ji oc an, the VEV of which is divergent due to 
the lack of the interaction term as stabilizer in Vn^afi) 
Eq. (jVI.gp . This reflects necessity of other operators such 
as the kinetic operator 9M('0V')t^^('0V') as counter terms. 
See Appendix IdI for further details. 

We depict our results mentioned above in Fig. |21 and 
Tables Ultall^: 



There could be a possible ambiguity of the estimate of the 
concrete number of Nf , Nc of the nontrivial window: It could 
arise from the ambiguity of the estimate of Kg'* in the im- 
proved ladder SD equation with different momentum identifi- 
cation for the running gauge coupling 10, 11|. Although, as we 
mentioned in Sect lllll our momentum identification Eq. 1111. 2i 
with nonlocal gauge fixing has much advantage over the sim- 



plest one |4j, 9f^_^.s-)DM ^ 9^4+i)D(™^^(^'' S^))' "^^^^ lan- 
dau gauge fixing, one might use the latter. Then one could 
get somewhat smaller values: [lOl Kg"' = 0.122, Kg"' = 0.146, 
and K^lo* ~ 0-163, compared with those in the present paper, 
Kgi' = D(D - 2)/[32{D - 1)]. This would shift the nontrivial 
(VF22) window for a fixed Nc to a somewhat smaller Nf region. For 
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(D = 6) 

TABLE I: Classification of models (D — 6). The entries la- 
beled by A, B and C stand for the models in the regions (A) 

KD > Kg", (B) K^T <; ^^ <; ^git ^nd (C) < I^D < I^d'^ , 

respectively. No entry denotes absence of the UVFP of the 
gauge coupling (?« . 



(A) The region kd > k™* - no finite continuum theory. 

As we discussed in Sec. IIIII the gauge coupling 
strength on the UVFP kjj is not a continuous pa- 
rameter, so that kd cannot be fine-tuned arbitrar- 
ily close to K^'* to make the dynamical mass finite 
in the continuum limit (A — > cx)). For example, 
models with D = 6,Nc = 3, Nf = 5, 6, 7 fall into 
this category. 



< KD < kS" 



nontrivial window. 



(B) The region k^'^ 

The dynamical mass can be made finite by fine- 
tuning of the four-fermion coupling. Once the dy- 
namical mass is made finite, the decay constant of 
the NG boson is also finite and so is the Yukawa 
coupling of the NG boson. The self-interaction 
of an and the Yukawa interaction remain nontriv- 
ial. We also note that the self-interaction of a^ 
is conformal symmetric. The renormalized four- 
fermion coupling gfj has a UVFP gji{oo) = g'^"* 
and the anomalous dimension 7^ is very large, 
i^ < 7m = J^(l + <-^) < h' + 1. For example, models 
with Z? = 6, iVc = 3, Nf ~ 4 fall into this category. 
The nontrivial window gets closed, k^'^ — > k^'*' 
when D —> 00. 

(C) The region < kd < n^'^ 
malizahle. 



trivial or nonrenor- 



Even if we made the dynamical mass finite by fine- 
tuning the four-fermion coupling, the decay con- 
stant does diverge in the continuum limit. The 



example, the model with D - 
nontrivial window. 



: 6, Nc = 3, N f = 3 would enter the 
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(D = 8) 

TABLE II: Classification of models (D = 8). The entries 
labeled by A, B and C stand for the models in the regions (A) 
Ko > k'W\ (B) k^"^ < kd < «;S" and (C) < ko < k^^, 
respectively. No entry denotes absence of the UVFP of the 
gauge coupling g, . For 2 < A^c < 5 the nontrivial window is 
not found. 



self-interaction of ur and the Yukawa interaction 
vanish and are trivial. The renormalized effective 
potential is not stabilized for gu > g"^"^. With- 
out introducing additional operators other than the 
four-fermion operator, the theory cannot be renor- 
malized. For example, models with D — 6, Nc = 
3,7V/ = 1,2,3 fall into this category. Note that 



,MT 



is reduced to the pure NJL point in the four- 



dimensional NJL model, 
D-> 4. 



^MT 



0, in the limit 



So far we focused our discussions on the theory 
with the gauge coupling exactly on the UVFP kd = 
C_f5*^nda- Each point on the critical fine g"^^ de- 
scribes a different theory (different Nc,Nf). Now that 
we have shown existence of the nontrivial window, we 
shall study the RG flows of (g'^dj,), gnifj.)) of a particu- 
lar theory specified by one point on the critical line. In 
order to obtain the RG fiows in the vicinity of the UVFP 
(5^(00) = gl, gu{oo) = g^"*^), we define variations Sg^ 
and Sgn as 



^5^(m) = 9^{t^)~ 9I, 
Sgnifi) ^ gnifi) - 5""- 

It is enough to study linearized RGEs [sl 
- (Sg^\ fSg^' 

ySgRj ygi 



(VI.24) 
(VI.25) 



dlnfi 



(VI.26) 
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(D = 10) 



since |Ai| - IA2I = 2(z/ - 1 - t/w) > for D > 6 (1/ > 2) 
in the nontrivial window {vco < 1). In order to estimate 
A, we need to determine the effective potential off the 
UVFP, g^ 7^ 5^, which can be obtained by solving the SD 
equation. Instead of doing this, here we simply replace 
the UVFP value gl in Eq. (|VI.4p with the running bulk 
gauge coupling g^: 



1 



1 



gnifi) ^{1+uj) 



:,^2 



fi = const. 



(VI.33) 



with 






(VI. 34) 



Differentiating Eq. IJVI.33|) with respect to In/i, we obtain 
a rough estimate of A as 



TABLE III: Classification of models {D = 10). The entries 
labeled by A, B and C stand for the models in the regions (A) 
KD > Kg", (B) K^'^ < KD < kS" and (C) < kd < k^'^ , 
respectively. No entry denotes absence of the UVFP of the 
gauge coupling g, . For Nc = 2 there is no UVFP. 



with 



M 



-2(i/-l) 
A -21^^ 



(VI.27) 



neglecting higher order terms of Sg^ and 6gB. The diago- 



nal components are obtained from Eqs. Ijll.5|) and (|VI.6|I . 
The off-diagonal one A is a constant to be determined 
later. The eigenvalues Xi {i — 1,2) and eigenvectors 
Vi , Wi are defined as 



vfM = XtvT, Mwi = \i 



W,; 



We obtain 



and 



Ai = -2(i^-l), A2 = -2i/w 



A 



Vl 



(VI.28) 



(VI.29) 



V2 = I I , (VI.30) 

2{iy-l- VLo) 



A = 2{v-l-vuj) 



l-uj g" 



9* 



(VI.35) 



The eigenvector wi is then parallel to the tangent of the 
critical line g"^^ . 

Once we choose the model, (iVc, iV/), the UVFP value 
(7* is fixed and so is the critical value of the four-fermion 
couphng g"^^. RG flows for a typical {Nc,Nf) are de- 
picted in Fig. O in the vicinity of the point Pqo = {g^ = 
From Eq. (jVI.32p we can see a gross fea- 



9*,9r 



9 



crit\ 



ture of RG flows independently of the value of A: The RG 
flows approach this point as increasing /i more rapidly in 
the direction of wi than that of W2 . There are four re- 
gions (I, II, III and IV) which are separated by two lines 
Li (parallel to wi) and L2 (parallel to W2). The D^SB 
takes place in the regions I and II. The line Li coincides 
with the tangent of the critical line. While the chiral 
symmetry is unbroken in the region III, the region IV 
is not precisely the symmetric phase: Although the NJL 
coupling in the region IV is not strong enough for the 
DxSB, the gauge coupling strength grows strong in the 
infrared there. We thus expect D^SB also in the re- 
gion IV, with a fermion dynamical mass being of order 

Al-^ (<C A) [nj which is a scale parameter of this gauge 
theory analogously to Aqcd.® Note that the theory is 
controlled by the UVFP Poo , where the composite fields 
a and -k enjoy nontrivial interactions. Thus the theory in 
the nontrivial window is consistently renormalized even 
off the UVFP 5 V 5*- 



Wl = 



-2{v- l-vuj) 



A 



W2 = 



We here note 



lAil > IA2 



\ ' ' 8 Thus Li as the border between region II and IV is not a phase 

boundary in the exact sense but has a similar feature in the 
sense that the dynamical mass grows rapidly from AL_: (instead 
of 0) to A as we cross the line from region IV to II. A similar 
phenomenon also takes place in the four-dimensional gauged NJL 

(VI. 32) model with the QCD-like gauge coupling. 
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FIG. 3: RG flows of a tfieory witfi D = Q,Nc 
Tlie point Poo denotes tlie UVFP, (5^(00) = 



= 3, iV/ = 4. 
ll,gR{oo) = 



rection toward tlie UV limit. The lines L\ and L2 correspond 
to the eigenvectors wi and W2, respectively. Li coincides with 
the tangent of the critical line g"^^ . 



theory defined at the critical point is only formal, since 
K£) cannot be fine-tuned arbitrarily close to the critical 
value [iO, [ij. The dynamical mass cannot be made 
finite in the continuum limit (A — > cxj). Then, even 
if we assume the existence of the UVFP beyond the 
truncated KK effective theory, the renormalizability of 
higher dimensional gauge theories with {Nc,Nf) such 
that Kb > K™* (region (A) in Fig. |2Jl is only formal ^. 
On the other hand, we do have a continuous parameter, 
the four-fermion coupling g, in the gauged NJL model 
and hence the nontrivial theory can be defined by fine- 
tuning g arbitrarily close to the critical line. Thus the 
inclusion of the four-fermion interactions may provide an 
interesting new possibility for higher dimensional gauge 
theories which were long considered nonrenormalizable 
and trivial theories based on the perturbation. 

Remarkably enough, we in fact found the nontrivial 
window (region (B) in Fig.|2) for {Nc,Nf) such that 






< Kd < K, 



(VII.2) 



VII. SUMMARY AND DISCUSSION 

We have discussed phase structure of the D- 
dimensional gauged NJL model with compactified extra 
(5(= D — A) dimensions on the TeV scale, based on the im- 
proved ladder SD equation with the running gauge cou- 
pling given by the truncated KK effective theory. Such 
a running (dimensionless) gauge coupling has a nontriv- 
ial UV fixed point (UVFP) and the theory behaves as a 
four-fermion theory coupled to a walking gauge theory 
with the gauge coupling almost constant near the value 
of the UVFP in a wide energy region above the compact- 
ification scale. The central assumption of this paper was 
that the existence of the UVFP is not an artifact of the 
KK effective theory and may have some reality in a more 
elaborate nonpcrturbative approach. 

Solving the SD equation in the bulk by setting the 
gauge coupling at the value of the UVFP as a reasonable 
approximation, we found the critical line similar to that 
of the four-dimensional gauged NJL model: 



D 

2 



1 



with 



D D-2 



l + Jl-nD/n-^ 



crit 



(VILl) 



which takes the same form as that 



^-D 32D-1' 

of the four-dimensional one with a fixed gauge coupling 
for D ^ 4 in the prefactor. In spite of the formal resem- 
blance to the four-dimensional case, however, the UVFP 
value Ku, given by Eq. HIII.6|I with Eqs. (|II.2|I and HII.7|I . 
is determined as a function of Nc and Nf and hence is not 
a free parameter but is a fixed quantity once we specify 
the model. 

Here we should emphasize the following: In the pure 
gauge theories with extra dimensions the "nontrivial" 






= 1- 



{DI2-1) 



(VIL3) 



where the four-fermion theory in the presence of walking 
gauge interactions becomes "nontrivial" and "renormal- 
izable" , similarly to the four-dimensional gauged NJL 
model, in the sense that the decay constant F.^ and 
couplings of composite tt and a as well as the induced 
Yukawa coupling of the fcrmion (with mass M) to those 
composites remain finite in the continuum limit K/M — *■ 
cxD. We explicitly performed renormalization of the ki- 
netic term of the composite bosons and the effective po- 
tential. The renormalized four-fermion coupling has a 
UVFP at gji{oo) = g"^^, where the theory has a large 
anomalous dimension 



In 



-?- 



1 + \/ 1 - kd/k 



.crit 
D 



(VII.4) 



The renormalized effective potential indeed has nontrivial 
interactions which are conformal invariant, thanks to the 
large anomalous dimension. 

It is rather surprising that even in higher dimensions 
the four-fermion operators become relevant thanks to the 
large anomalous dimension, D/2 >_-fm > D/2 — 1. In 
fact, the dynamical dimension of (i/'V')^ operator is rel- 
evant, dim((V'V')^) < D. On the other hand, the {'ipip)* 
operator and the kinetic term dM{'4''4})d^ {'ipi') are ir- 
relevant, ^{^{{Tpi^f) =2D- 2{D - 2)^1-kz5/k™* > 
D, and dim(9M(V'V')9^(V'V')) = D + 2 - {D - 



^ Nevertheless, such a theory can still be useful for model building 
as an effective theory with finite cutoff A. See, for example, the 
tMAC analysis of the D = 8 Top Mode Standard Model with 
extra dimension |12| . 
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(D = 4) 



FIG. 4: RG flows of a theory for D = A. The arrows of the 
flows indicate the direction toward the UV limit. The line Li 
corresponds to the eigenvector wi (tangent of the critical line 
of the four- dimensional gauged NJL model with flxed gauge 
coupling |2a.l27l|'). The black circle is the pure NJL point. 



TABLE IV: Classification of models (D = 4). The entries 
labeled by B and C stand for the models in the regions (B) 
A = 18Cf/(11A^c - 2Nf) > 1 and (G) ^ < 1, respectively. 



asymptotic behavior of the mass function |2E 



2)-\/l — k_d/k™' > D. Thus the theory can be renor- 
malized without operators other than (ipip)'^. We further 
gave the RG flow off the UVFP in Fig. O which is consis- 
tent with our renormahzation performed on the UVFP. 

On the other hand, for {Nc, Nf) such that < kjj < 






(region (C) in Fig. [SJ, the decay constant F^r di- 
verges even if we fine-tune the four-fermion couphng close 
to the critical value g"^^ to make finite the dynamical 
mass of the fermion. At the MT point kd = k^'^ the 



divergence of F^ is logarithmic, F^ ~ In A, which is sim- 
ilar to the situation at the pure NJL point of the four- 
dimensional case. The MT point is reduced to the pure 
NJL point K^^ — > for _D ^ 4. Without introducing 
other operators such as {"ipipy, 9m(^^)5''^(V'V'), etc., the 
renormahzation breaks down in the region (C) 

At this point it is worthwhile mentioning the formal 
limit £) — » 4 of the nontrivial window. It was shown p35l 
ISo . 1371 IsM l39| that the four-dimensional gauged NJL 
model becomes renormalizable and nontrivial, when the 
one-loop gauge coupling runs as specified by: 



A 



6Cf 



> 1, 



(VII.5) 



Cf 



N?. 



2Nr_ 



-b = 



llNr-2N 



where A measures the running speed of the coupling [22, 
|33|:A ::^ 1 is the walking theory and A — > cx3 (6 — > 0) cor- 
responds to the non-running (standing limit of walking 
coupling) case with vanishing beta function. The condi- 
tion Eq. IjVILSp was first obtained |35j by evaluating the 
F^, taking account of the logarithmic corrections to the 



T,{x) r^ ciEo 



'"'i 



(VII.6) 



instead of Eq. IJIII.36P , and then in the RG group anal- 
ysis in the language of the equivalent gauged Yukawa 
model "SSISTL l3a | and also in the nonperturbative RG 
eauation|i39j. 

Here we show that the formal limit of I? ^ 4 of the 
nontrivial window in this paper coincides with that char- 
acterized by precisely the same condition as Eq. HVII. 5|l , 
now without logarithmic corrections which were rather 
delicate factors in the previous arguments in the four- 
dimensional model. Instead we consider = 4 + 6 with 
< (5 <C 1 in the spirit of e expansion [SJ, |5a|. The 
nontrivial window is defined by 



VU! < 1, 



which reads 



^D 



>S. 



(VII.7) 



(VII.^ 



On the other hand, from Eqs. (|II.7|I and ljIII.6|) we have 



Kp 

t^crit 



6Cf 
~6" 



AS 



(VII.9) 



Then we obtain Eq. IIVII.51) . A > 1. The nontrivial win- 
dow in the formal limit D ^ 4 {6 ^ 0) thus coincides 
with the condition of the renormalizability/nontriviality 
of the four-dimensional gauged NJL models, namely 
those coupled to the (moderately) walking gauge the- 
ory with A > 1. Combining the condition of asymptotic 
freedom Nf < llNc/2 and A > 1, we thus have theories 
within the nontrivial window in D —f 4: 



9 1 

27Ve 



<N 



f 



Nc < -N,. 



(VII. 10) 
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See Table |^ for the theories satisfying Eq. HVII.10|) (or 
Eq. (|V11.5|l ). On the other hand, the trivial region (C) 
{0 < kd < K^'^) coincides with the condition of the 
triviality/nonrcnornializability, A < 1, in the D —* A 
limit. Thus in the Z? — + 4 limit where k^"^ — *■ and 
(for A < oo)!'', the UVFP's for both regions 



KD 



^MT 



< KD < K™*) 



(C) {0 < KD < K^^) and (B) (k^ 
the phase diagram of Fig. |21 shrink to a single point 
of the pure NJL point and hence the distinction among 
them is not obviously visible compared with the case 
with extra dimensions. In other words, the renormaliz- 
ability/nontriviality of the four-dimensional gauged NJL 
model with A> 1 gauge theories [sl IMIstL l3& .SQJ is a 
four-dimensional manifestation of the nontrivial window, 
in sharp contrast to the theory with A < 1. Thus the 
nontrivial window is not a peculiarity nor an artifact of 
the extra dimensions but is rather a universal feature of 
the gauged NJL model. 

The RG flows of the nontrivial window in this limit 
are only in regions II and IV of Fig. |21 and are similar to 
Fig. but with a crucial difference: Since |Ai| < IA2I in 
this limit, the flows first approach Li instead of L2 and 
then converge toward the pure NJL point (Fig.0}. This 
is actually consistent with t he g ross feature obtained in 
the four-dimensional model. JS^, |33, |53 

In conclusion, we have shown that the four-fermion 
theories in Z) > 4 dimensions are renormalizable and 
nontrivial, when coupled to D-dimensional walking gauge 
theories, with the parameters {Nc, Nf) being in the non- 
trivial window, which we assumed have a nontrivial ultra- 
violet fixed point as given by the truncated KK effective 
theory. The fact that the four-fermion theory in dimen- 
sions D < 4, lower than 4, is renormalizable and nontriv- 
ial has been known for long time since Wilson |5a| (For 
a proof on the D = 3 theory, see Ref.|53) and was also 
known for the arbitrary £)(< 4) dimensions 1371 l59l l60l| . 
Further in £> = 4 it was known 111131131138^39] that 
the four-fermion theory when coupled with certain walk- 
ing gauge theories is renormalizable and nontrivial. Here 
we have shown for the first time that the four-fermion 
theory in I? > 4, higher than 4-, also shares the same 
feature when coupled to walking gauge theories charac- 
terized by the nontrivial window. 

This renormalizability/nontriviality is not just a for- 
mal matter, but implies cutoff ins ensitivity oi the physics 
prediction, which may be useful for model buildings even 
with a finite cutoff. This is rather surprising, since the 
theory is dominated by the dynamics in the ultravio- 
let region where both the gauge and the four-fermion 
interactions become strongly coupled. Explicit model 
buildings based on this observation were in fact at- 



^^ If one formally changes the theory [Nc, Nf) depending upon the 
dimension <5(— > 0) in such a way that A ^ 4A/(5 — > 00, one could 
have a limit arriving at the four-dimensional gauged NJL model 
with fixed (standing) gauge coupling A: Ku/reg'' — > 4A (0 < A < 
A<;"t = 1/4), 



tempted for the four-dimensional case ( "Top Mode Walk- 
ing GUTs") 13^,61,]. 

The phase structure of the gauged NJL model in the 
bulk with compactified extra dimensions may be use- 
ful for model buildings such as the Top Mode Standard 
Model, bulk technicolor, etc.. For instance, the D = 6 
Standard Model gauge interactions are not enough to 
trigger the top quark condensate [ll [Uj 12], while the 
gauged NJL model can do work due to the additional 
attractive interactions as in the case of the original Top 
Mode Standard Model. Although introduction of such 
an ad hoc four-fermion interaction may be less attractive 
than the scenario without it, the KK modes of the top 
quark still naturally reduce the top quark mass predic- 
tion of the original model as was emphasized in Ref. J3] ■ 
Since the electroweak symmetry breaking is still a cen- 
tral mystery of the modern particle theory, it would be 
useful to consider all possible dynamical scenarios before 
the LHG will take off. The gauged NJL model with extra 
dimensions may become one of the dormant volcano to 
give an explosion into rich TcV physics. 
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APPENDIX A: PROPAGATOR OF THE 
COMPOSITE SCALAR 

We have calculated the auxiliary field propagator in 
Sec. lIIIDl bv using the resummation technique 48]. We 
may take yet another choice proposed in Ref. [43 in 
which we obtain the scalar propagator as the lowest or- 
der of the Ghebyshev expansion. While the resummation 
technique is operative in a weak gauge coupling region, 
the Ghebyshev expansion method is valid in the whole 
region. In four dimensions the results of the Gheby- 
shev expansion method are similar to those of the re- 
summation technique. [43 We here demonstrate that the 
wave function renormalization constant for the auxiliary 
field agrees with the result of the resummation technique 
shown in Eq. (|III.67ll even in extra dimensions. 

The composite scalar propagator is written in terms of 
the vacuum polarization function lis for the composite 
scalar: 



^D-\q) = -^+^s{q), 



(A.l) 



where the vacuum polarization function is given by 

ns(g) - N,Nf J -^^ tr [S{k + q)xs{k + q, k)S{k)] , 

(A.2) 



with xs being the Bethe-Salpeter (BS) amphtude for the 
scalar. Neglecting the explicit breaking of the D dimen- 
sional Lorentz symmetry owing to the compactification, 
we can decompose generally the BS amplitude xs into 
four Lorentz scalar functions Fi {i — 1,2,3,4), 



xs{p + q.p) = 

Fi{p + q,p) + {ii> - i>i) F2{p + q,p) 
+ {^ + ^)F3ip + q,p)+'^Fi{p + q,p), 



(A.3) 
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scalar bound state, respectively. Although one might sus- 
pect that more structure functions are needed in D{> 4) 
dimensions, other terms are zero or reduced into the 
above four, because the BS amplitude includes only two 
momenta. We analyze the scalar propagator and the BS 
equation in the symmetric phase, i.e.. 



S{p) = Sq{p) 



(A.4) 



where p and q denote momenta of the fermion and the We then obtain the vacuum polarization function as 



-,. i,r 1 /^ I d^k -k-(k + q) , , [ dPk ik-qf-k'^q^ , 



i{2Tr)D B{k + qf '^ < H,"-; < -J -^^^^d p(fc + g)2 
On the other hand, the BS equation under the ladder approximation is 

d^fc 



(A.5) 



xsip + q,p) 



Nt2 



(2^) 



D I 



iT^r''']Sik + q)xsik + g, k)S{k)[-iT''r'']g(^^s)DDMN{p - k), (A.6) 



where 1 is the identity matrix and Dmn is the propagator of the gauge boson. In the symmetric phase the BS 
equation is obviously decomposed into two simultaneous equations for i^i.2 and for -F3.4. Since the structure functions 
F^ and F4 do not contribute to lis as shown in Eq. (jA.5p . we neglect the two hereafter. The BS equations for Fi and 
F2 are given by 



F^{p + q,p) = 1 + {D-1 + S,)Cf91 



F2ip + q,p) 



d^k 



-k- (k + q) 



i{2Tr)D fc2 {k + qf [-[p - A:)2]^/2 
d^k {k ■ qf - fc2q2 



— Fi(fc + (7,fc) 



Cpgl f d^k S(p,g,fc) 

2[{p-qy-p^q^]Ji{2T:)Dk^k + q)'^[-{p~k)^]D/2~i ' 

, Cpgl f d^k -k-ik + q)E{p,q,k) 

{p-q^-p'^q^J i(27r)^fc2(fc + g)2[_(p_A:)2]W2-i 2 1 +9, ) 



with 



E{p,q,k) = {D - 3 + 0[q\p ■ k) - {p ■ q){q ■ k)] 



2(1-0 
-(p-fc)2 



{-q^[k ■ (p - k)] [p-{p~ k)] - [q ■ {k-p)] [p\q ■ k) - k^p ■ q)]} 



r 



(A.7) 



(A., 



(A.9) 



where the bulk gauge coupling was replaced as in depend on three angles a, (3, 7 defined by 

Pe ■ qs 
cos a — 



Eq. Ijlll.211 with g^ = g^. The gauge fixing parameter 
^ should be taken as 

{D-l){D-i) 



\PE\\qE\ 



r. PE-kE qE ■ kE 

cos p = - — — — r , cos 7 : 



\PE\\kE\'' 



D 



for consistency with the ladder approximation. 

Let us expand Fi, F2, and the integral kernels of lis 
and the BS equations by series of Gegenbauer polyno- 
mials, which is generalization of the Chebyshev polyno- 
mials. At this stage the BS equations HA.7|I and ljA.8|) 



qE\\kE\ ' 
(A.ll) 
as well as an infinite chain of harmonics /„ and g„ for ex- 
pansion of the Gegenbauer polynomials C^ (cos a) , where 



Fi{p + q,p) = Y,Upl,ql)Cl{cosa) (A.12) 



n=0 
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and 



F2{p + q,p) = y^Qn{p%,q%)Cl'^{cosa), 



n=0 



(A.13) 



so that it is quite difficult to solve analytically the equa- 
tions. We here note that only the BS equation for the 
harmonics /o contains an inhomogeneous term (the con- 
stant unity), because the Gegenbauer polynomials satisfy 



dd sin^'' e C^ (cos d) C,"; (cos 9) = w,, S„ 



with 



and 



w„ 



7rr{n + 2v) 



^{n + iy)nl[r{v)Y 



C^{cos0)^l. 



(A.14) 



(A.15) 



(A.16) 



The inhomogeneous term certainly gives main contribu- 
tions to the BS equations. Hence we assume that the BS 
equation for /q is well approximated by the closed form 
of /o, neglecting the effects of /„(« > 1) and (?„(« > 0). 
The harmonics /«(« > 1) and gn{n > 0) are iteratively 
determined after /q is computed. We then obtain the BS 
equation for /o as 

4( D — l) f^^ 
/o(s,t) = 1 + -^^-^^ — -i^D / duu^-'^ Ko{s,t,u)fo{u,t), 



where s = p'^ , t = q^, u — k'i,, and 



^B' 



Ko{s, t, u) = an{u, t)bo{s, u) 



(A.17) 
(A.18) 



with 

ao{u,t) 
1 



d7sm^''7C^(cos7)— - — ■ --,(A.19) 

Wo Jo [kE + qEr 



6o(s,m) 



C^(cos/3) 






(A.20) 



We can easily perform the angular integrations and find 
1 



ao(w,i) 



bo{s,u) 



u — t ^( D D min(u,t) 



max(M, t) 



2 2 max(u, t) 



[max(s,u) 



(A.21) 

(A.22) 



The vacuum polarization function is written in terms of 
/o as 

ns(i) = 2*rjNDAiVciV/ / duu^-^ao{u,t)h{u,t). 
Jo 

(A.23) 



Noting that the integral of R.H.S. in the zeroth order BS 
equation IjA.lTII at s = A^ is just same as in Eq. (|A.23|I . 
we can rewrite the composite scalar propagator without 
the integral as follows: 



[Fuv(A2,i)-l],(A.24) 



1 
9 



4(D-1) 



where we divided the zeroth order harmonics /o into two 
parts of the IR and UV regions, 

/o(s,i) = Fm{s,t)e{t -s) + Fuv(s,t)0(s - t). (A.25) 

The BS amplitude in zero momentum transfer, i.e., the 
effective yukawa coupling Fs, is also obtained through 
-Fuv, 



^s{Pe) = XsiP,p) - FuY{p\,q\ = 0). 



(A.26) 



Let us solve the zeroth order BS equation (|A.17|) which 
is equivalent to a set of the second order differential equa- 
tion, 



^/o(s,t) + — — /o(s,i) + ^ao{s,t)fo{s,t) - 0, 

(A.27) 
and the IRBC 



9s2 



s^/^^/o(^,t) 



= 0, 



s^O 



and the UVBC 

s^/o(s,i) + (^/2-l)(/o(s,t)-l) 



with 



(A.28) 



= 0, (A.29) 



s^A2 



_ 2{D-l){D-2) {D~2f{l-u^) 
kd — ;:; K-D ~ —. • (,A.3L)j 



D 



16 



We now replace the hypergeometric function in ao{u,t) 
by unity, i.e., 



ao{u,t) 



u — t 
max(M, t) 



+ 1 



(A.31) 



We confirm later on that the approximation works well 
in D ~ 6. Within the approximation the differential 
equations for the UV and IR parts are 

d iy + 1 d ko f t\ 



and 

^FiR + ^-FiR-t-— FlR-0, 

respectively, where we used 

v^D/2- 1. 



(A.32) 
(A.33) 

(A.34) 
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We can solve analytically Eqs. IJA.32|) and ljA.33|) and find 



FvYis,t) = ( ^ 



2k£,- 



+C2liyuj \ r^KD 



t 



(A.35) 



for the UV part, and 
Fm{s.t) = 



c:iJ^ W2k_d 



t 



+C4Y1, [ \j2kD- 



(A.36) 



for the IR part, where J^, Y^ and I±vi^ represent the 
Bessel functions of first and second kind, and modified 
one, respectively. For uj — 1/v we should use the modified 
Bessel function of the third kind, Ki, instead of /_i. The 
IRBC (TOl leads to 



C4 = 0. 



(A.37) 



Other coefficients ci,C2, and C3 are obtained from the 
UVBC (|A.29|) and smoothness conditions for /o(s,t) at 
s = t. 



c(^-- 



-(- 

2v I A2 



Viv 



2kjjt 






2k£)t 

~A^ 



(A.45) 



The prime X'{z) denotes dX/dz. We also note the rela- 
tion 



/;(z)/_,(z) - Uz)lL,{z) - ^^^I^M, (A.46) 



The BS amplitude is thus given by 



^ih(M)=.-M^--)G)'^4V^MA.47) 



Fvw{s,t) 



2 sin(7ri^w) "^ \A'^''^J \s 



t 



1 {^) I-i- 



2kut 



-7 (-w) I^ 



"Ihjjt 



(A.48) 



Fir(s -^ t,t) = -Fuv(s —> t,'t)^ (A. 38) The effective yukawa coupling Fg is also obtained as 



and 

d_ 
'd's 

as follows: 



Fm{s.t) 



ds 



FvYis,t) 



Ci = Ci 1 ^,w 

7r7 (oj) 



2sm{Tn'uj) \A2 



C2 = ci 1 ^,-^ 



C3 ^ V 



where we defined 



A^'^ 



A2 



2 sin(7ri/cj) 



7MC( ^2^ 



(A.39) 



(A.40) 

(A.41) 
(A.42) 



--fi-^)C\ ^,w 



, (A.43) 



7 (w) EE V2^ [j>^ (\/2^) /L (\/2^) 



JM 72^) /,^ ( V2^) , (A.44) 



r.bl)=Fuv(p|,g| = o) = ^('J') ; (A.49) 



which agrees with Eq. Ijin.53(l calculated through the re- 
summation technique. 

In the case of I? = 6, we can solve the BS equation 
without using the approximation for cq represented in 
Eq. HA.31|I . The BS equations for Fir and Fyy are 






d ^ , s 3 d ^ , 

— Fuy{.S,t) + -—F^y{s,t) 



ds 



s ds 



respectively. Noting that Eq. IJA.51|) is rewritten as 



d Id f I 2t e\ 



ds 



s ds 



3s 6 



(A.52) 
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with s = 1/s, we obtain the solutions 

Fuv(s,t) = 

fket 




+C4 exp 



where i-Fi (01,7; z) and ^>{a,"f;z) denote the confluent 
hypergeometric function and its $ function, respectively, 
We here defined 




7± = 1 T 2^/1-^6, (A.55) 

aiR=--2W-f. A.56 

2 V 



We can confirm that Fuv is real, i.e., F^ 



uv 



F, 



uv, 



by using the Kummer's transformation iFi [a,j; z) — 
e\Fi (7 - a, 7; -z). The IRBC (TOS^ leads to C4 = 
owing to the behavior of the confluent hypergeometric ^ 
function around zero 

*(a,l + n;z)~ Ilizill^-" 



(-1)'^ 



n\ T{a — n) 



r(«) 

lFi (a, n + 1; z) In z, (z - 0). (A.57) 



Other coefficients Ci , C2 and C3 are obtained by us- 
ing the UVBC (JA.29P and continuous conditions HA.38|I - 
(jA.39p . We show the difference between the solutions 
with or without the approximation ljA.31|) . (See Figs.[Sl 
and 0) The approximation works well within a few 
%. We expect that the approximation is also valid in 
£1 = 8,10, •••. 

We now determine the wave function renormalization 
constant for the auxiliary field a. We expand the function 
Fuv(A^, g|) in Eq. HA.48|) over ^l/A^. Depending on the 
values of vlo, three expressions are obtained: 



iD-\ql) = -2^ NcNf^^YiA 
'I 1 



xA 



D-2 



9 9' 



crit 



Ai 



(A., 



for < w < 1/i^, and 

iD-\ql) = -2*7Vc7V/r!NDA 
xA^-2 



1 

9 


1 nl^ /A2\] 



(A.59) 



0.005 




FIG. 5: Error of the zeroth order BS amplitude /o arising from 
the approximation <A.3H . The curves show {f^^^ — /o)//o 
at s = A^ in the unit of A — 1, where /o and f^^^ denote 
the solutions in Eqs. ljA.53^ - l|A.54^ without the approxima- 



and the approximate expressions in Eqs. 



tion 

1A.48I I. respectively. The dot-dashed, solid, and dashed lines 
represent graphs for Kg = we/Ke"* = 0.02,20/33,0.98, re- 
spectively. (The ACDH model in D = 6, iV^ = 3, Nf = 2, 
corresponds to k^ = 20/33.) 




FIG. 6: Error of the zeroth order BS amplitude /o arising 
from the approximation 1A.31II . The curved surface shows 
ifo''^ - /o)//o for fte = 20/33 {N^ = 3,Nf = 2) in the unit o f 
A = 1, where /o and /g^'' denote the solutions in Eqs. HA. 5311 - 
IA.54I1 without the approximation IIA.311 1 and the approxi- 

respectively. 



mate expressions in Eqs. 



for uj = \/v, and 

zD;l(4)--2*iVeiV/f2NDA 

■f 1 



xA 



D-2 



9 9 



A ^^ 



(A.60) 
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for l/v < oj < 1, where we defined 

_ j{—Lo) r(l — i/uj) 16a; 



Ai 



7(0;) Til + iyuj)iy{l+uj)3{l-Lu) \ 2 



kd 



A2 = 



2v' 



[v + iy 



and 



A. 



{1 + ujY{1-vuj)' 



is unciianged even in the broken phase, as long as neglect- 
ing effects of the mass term. We thus obtain the wave 
function renormalization constant for a as 



(A.61) 


f 




2i/(l-(.j) 




for < cj < - 




z-'- 


) 






(A.62) 

< 


z-' = 


/A^ 


2(^-1) 


In 




for tj = — 


(A.63) 


Z-' - 


/A^ 


2(^-1) 






for - < w < 1 

V 



(A.65) 



For w = 1/;/ we also used the behavior of the Bessel ^\\\i 
function Ki{z) around z ~ 0, 

Xi(z)- - + |lnz, z-0. (A.64) 



an - Z-i/V. 



(A.66) 



Although we have calculated in the symmetric phase so The results are consistent with the analysis of the resum- 
far, we expect that the momentum dependence of D^^{q) mation technique in Sec. IIII Dl 
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APPENDIX B: LINEARIZED SD EQUATION 

In this appendix, we summarize our results on the linearized SD equation including the cases with kd > k^'' and 
K£) = K^'^ . We obtain the asymptotic behavior of the mass function as follows: 



E(^ 
So 



X \ ^ { \ X 

2|ci| ( ^ j sin(-iyJ)ln^ 



T{D/2) 



r(i./2)r(i + i./2) Vs2 



X 
^0 



dl 



^0 



'-4<'''-"t'°52'' '- 



-l(l-c^) 



'0. 



for 



for 



for 



for 



for 



KD > Kg'* 



Kd — Kjj , 



Kjj < Kd < K^ , 






0< KD < K^^, 



(B.l) 



Kd = 0, 



where 



Cl 



r{D/2)r{wCj) 



r(j/(i + ij))/2)r(i + j/(i + ic:;)/2) ^ 



dl = c*, 



-ci/di 



Kd 



1, 



(B.2) 



Cl 



r(L>/2)r(i/w) 



r(zy(l+a;)/2)r(l + j.(l+a;)/2) 



, dl 



r(i:)/2)r(-j/a;) 



r(t/(l - w)/2)r(l + i^(l - a;)/2) 



' ^2 = --77 ryci, B.3 



and 



ho = -2j-2yj{iy/2)-^, VW = ^lnr(z). 



(B.4) 
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In Eq. HB.4|) . 7 is the Euler's constant, 7 ~ 0.5772. 

Substituting Eq. (jB.ip for the UVBC (|III.18|I with toq = 0, we obtain the gap equation 



sin ( vuj In \- 

So 



\n-^ = h 



16 



V g 



4wdi 



So 

A 



2i/ci; 



1-Cj2 



Cl 



1+ci; 



dl 



So 

A 



2f cj " 



^ — = 1 _ t.5"" ( 1 - 



4C2 



(^-1)M^I]2^^^A2 



4i. 5; A2 E2- 



A2 



^(^--^)c,+(i-^c.i^ 



= 1 - 



L 9 



- 1 A2 ^ 



for 



for 



for 



KD > Kg'* 



KjJ Kjy , 



for Kjy < K£, < K 



D 1 



for KD ~ Kd I 



for 0<kd<k'S, 



kd = 0. 



(B.5) 



with 



tan 6*' = 4^ ^^^ r^'^- 



(B.6) 



The scahng relation is then found as 



So 

A 



exp 



exp 



exp( -) , 



hn 1 



exp 



= < 



Cl 1-uj- ( ^ ff 



'di iiiJ 



4 


1 


\ ^ 


9 


crit \ ■ 


1 

2f 


9 )\ 





Cc--- 1- 



for 



for 



for 



Kd > Hi 



K]J Kjy , 



K£) < K]J < K£) , 



for 0<kd<kd , 



(B.7) 



where we used g ~ g'^" • For kd = kj^ the scahng relation cannot be written by elementary functions: 



Sg, A^_ 4 g-^ 
A2 Eg (1 + ^)2 (.^ 



(B., 
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The chiral condensation is calculated from Eq. (jIII.lSp : 



|ci I V 1 + w^ ( — J sin ivLjln— — \- 9 + arctanw 



r{D/2) /Soy /A ho 1 

T{iy/2)Til + u/2) \ A J \ ^q 2 ly 






for 



for 



KD > Kg'* 



KU — Kjj , 






MT 



for k'^' <kd< Hd, 



(l + iy) /S 



-ci 



iz-l 



2i/ V A ^ 

1, 



1 y2 a2 



for 






2 -^^A 



1+LU 1 



2 V \ K 



- C2 



i^(l-c^)+2 



for < KD < K^"^, 

for K]j — 0. 



(B.9) 



Eq. HB.9|I reads 



2fT 



2(i^ + l) a 



2 2w , . / 2c7^ a 
-a+z diA' ^ 



l+T 



E(A2) = < 



1+W l+LU \l+wA 

2iy i^-l , [2vc:^ oY""^ ^ A 
a CiA ^—-r\ In—, 



2 2 / 2c{' a 



for 



KD = kS'\ 



for K^T < K^ < ^crit^ 



for 






for < KD < Kd^, 

for kd = 0. 



(B.IO) 



For ky) > K^'* we do not obtain a handy formula. 



APPENDIX C: POWER EXPANSION SOLUTION 



We summarize results on the power expansion solution for ko < k^J'*. The asymptotic behavior of the mass function 
is given by 



S(.t) 



-ii/(i-(^) 



^1 I ^ 



-il/(l + Uj) 



for K^^ < KD < K'j^'^ 



= < 






■^0 



Cl> J.2 



(C.l) 



■^0 



-il/(l-(^) 



with 



C2 , 2 

^1 



C2 



-|i/(l-w)-l 



for < K_D < K™, 



Z/2(l_c^2) 



4(j/ + 1 - iyLu){2iyuj -V -I) 






(C.2) 



where ci and di cannot be determined in the power expansion method. 

Substituting Eq. (|C.1() for the UVBC Hill. 18(1 with mp = 0, we obtain the gap equation 
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5"''* _ , Audi 



So 

A 



1-^ /Sn 



2uijj ' 






9 



v + l 



5"''^ _ 4(i/ + 1 - vuj)c2 



^crit 

9 



So 

A 



2;y(l-w)+2 



z/(l-a;2) 



Cl 



i/(3ccj~l)-2 



..(l+o;) "Ma 



C2 



So 



2y(l-ui)+2 ■ 



for K^^ < K_D < K™*) 



for 






for Q < KD < K^- 



(C.3) 



The scahng relation is then found as 



So 
A 



2z^a; — j^ — 1 



1- 



1- 



for n^j^ < KD < kS , 



for < K_D < K^ , 



(C.4) 



where we used g ~ 5*^" . For no = Hjj the scahng relation cannot be written by elementary functions: 



Eo\''+\ A2 Z.+ 1 g~^-2l^±^ 



S2 2l^2c2„crit _ (^-1)^ 

i» y 16!^^ 



(C.5) 



The chiral condensation is calculated from Eq. pil.l 



(1+c^) /So 



iy(l-(.j) 






^ _ I (3Z/+1) /So 

T^ — \ Cl ^ 



4. ^Ma 



2^Q,, , 1^ /v„\ ^^ A2 



{u-lY{iv+l) 3 /So 
32iy(iy+l) "^n A 



c?(^) In-^, for 






y2' 



crit 



for K^^ < KD < Kd^ 



kd — Hjj , 



for Q < KD < k^d'- 



(C.6) 



Eq. l(Cl6)l reads 



S(A") = 



2 2cj , . /2c7^ cr\^^ 
-(T + ^ diA ' ^ ' 



= < 



l + uj \+uj V1 + cjA 



Av V -I ^ f Av a\" , K 
-a A In- 



3i^+l v + 1 \'iv+lK a 



v{\-uj) 



■A 



2 a 



I 



for K]3 < kd < n' 



D J 



for 



for 



I^D — f^D : 



< KD < K^- 



(C.7) 



These formulae are useful for calculation of the effective potential. 
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APPENDIX D: MORE ON 
NONRENORMALIZABLE REGION 

In Sec. I VII we studied the renormalization of the effec- 
tive potential, keeping the kinetic term of the auxihary 
field finite. It should be so in the viewpoint of renor- 
malization of the effective action. For < k_d < K^f^, 
however, the theory cannot get into the "broken phase" 
gn > g"^^, because the "renormalized" potential ljVI.9|) 
becomes unbounded from below. In this appendix, we 
attempt another procedure in which the VEV of ct/j re- 
mains finite independently of the value of kjj. For such a 
purpose we may admit to incorporate any required oper- 
ators. Can be found a necessary and sufficient operator 
set? Before details, we jump to conclusions: It depends 
on approximation methods of the ladder SD equation. 

The dynamical mass should be finite, so that 
Eq. IJIII.15(I yields the wave function renormalization con- 
stant of a as 






2i/(l-w) 



iO<KD< Kg'*) 



(D.l) 



for the whole region. (Compare Eq. l|D.Hl with 
Eq. (jIII.67|l .) By using Eq. ID.ip . we obtain the renor- 
malized Yukawa vertex 



rf)(-,^) 



2\ --1(1-") 

^ , (0 < ACi, < <*), 



/^ 



(D.2) 



in the continuum limit A ^ cx). (See also Eq. (|111.69|l .) 
On the other hand, the "renormalized" propagator of a 
behaves as 



iD-Up)-^D-UO)'^ 



<y(R) 



^(i?,)V 



A 



2(iJU-l) 



for 0< kd < nW and 



^D~lR)ip)-^D-lR)i^)^^?^ (^ 



M" l^^ ] (D.3) 



^y{^^ (°-^) 



for kd = f^"^ ■ Thus the kinetic term of a diverges in 
the continuum limit A ^ oo for < k^ < k^^ ■ We 
need to add the operator dM{^'4')d^ {'ipi^) in the effective 
action in order to absorb the divergence. Moreover, the 
decay constant is also divergent in the region, so that the 
interactions of the NG boson vanishes in A — s- oo. 

We investigate the effective potential for < kd < 
K^"^ separately in two approximations. 



1. Linearized approximation 

The scaling relation for < kd < KjP' is given by 



Eq ex Aw 1 



5' 



crit 



9 



(D.5) 



(See Eq. IjIV.iyp .') The four-fermion coupling g is thus 
renormalized as 



1 



1 

^crit 



A^ = 



1 



1 



Ai 



(D.6) 



5(A) 5™V'' V5i?(M) r 

so as to make the dynamical mass Eq finite. This is the 
very same as Eq. IJVI.11|) . 

The bare effective potential has been already obtained 
in Eq. (|V.16(I . We rewrite Eq. (|V.16() in terms of ajj and 

ImoaR f A^ 



Vian) ex ^^ 



+A, 



9 Ai 
2 \9R 



1 



M 



2+7 



4 



2(iyw-l) 



2(iya)-l) 



(D.7) 



We renormalize the bare mass term, i.e., the first line of 
Eq. HD.7|I as follows: 



5 (A) 9r{i^) 



(D., 



which is the same as Eq. (jVI.3|) . The beta function and 
anomalous dimensions are then found. 



f3{9R) = 2gfl 1 



9R 



(D.9) 



and 



7m (ff/ 



P{9r) 
9R 



+ v{l+uj), 



v{l + uj) - 2 + 2 



9R 



9^ 



(D.IO) 



respectively. At the UVFP of gn the anomalous dimen- 
sion is 7m = v{l + Lj). Therefore the mismatch of the 
A dependence between toq and {{'>p4')o) is resolved. (See 
also Eq. fTTm .) 

Troublesome are the aj^ term and interaction term {A2 
term) in Eq. IjD.Tp . The divergence still remains after the 
"renormalization" of the four-fermion term. In addition, 
the operator having the fractional power of (ip'tp) is re- 
quired. 



2. Power expansion method 



The scaling relation is given by 

1 

„crit \ 2u^ 



A 1- 



9 



En~ <^ 



9 



Crit \ 2iy{l~ui} + 2 



{Kjj < kd < k 



D ^ 



PEn 



A ( 1 - 2 — ) ^ (0 < KD < K^' 



(D.ll) 
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(See Eq. (|IV. 33^ .1 In order to keep Sq finite, we renor- 
mafize the gap equation, 



for K,^ < Ko < K™*, and 



1 



A 



2i^aj 



g{A) 5"it 
for K^ < kd < K^J'* and 

1 1 

ff(Ay~5^ 
1 



9Ri^^) 5"'* 



A 2i/(l-w)+2 



,2i/w 



(D.12) 



gRifj-) 



M 



2u(l-Lj)+2 



(D.13) 



for < Kd < K^. By rewriting the bare effective po- 
tential in Eqs. (|V.20(I and (|V.21(I in terms of an. and gn, 
we obtain 



V{aR) ex fj. 



D 



1 rnpaji /A 
5 M^ VM 



1^(1+^^) 



-(-- — 

2 V.gj^ 5"'* 



1 A 4 



/i- 



+Ai 



CTi? 



F(afl,) « A^^ 


1 moan 
9 M^ 


V/v 




1 / 1 


1 \ al /Ay(2..>-.-i) 


^2 Ufl 


5"'V ^l^\^i) 




4 /^x 2(27.^-1.-1) 

) w 


: 



(D.15) 



for < Kd < K^ . The bare mass term is then renor- 
mafized 



"^O(^) At'd+c^) ^ '71fl(^) ^(1+^) 



ff(A) 



5ii(M) 



(D.16) 



for < Kd < Kn'*- The beta function and anomalous 



(D.14) dimensions are 
I 



and 



7m (5i?) 



P{9B.) 



9B. 



'^"^an. ( 1 - ^ 



2[i.(l-c^) + l]5i^ 1 



+ iy{l+uj) = 



9R 



\-uj + 2uj 



9r 



(k£^ <kd< «B'*) 



[0<KD< kPJ') 



3uLu ~iy -2 + 2[i/(l - w) + 1] 



9R 



PE 



< Kd < k' 



D 



{0<KD< kI^) 



(D.17) 



(D.18) 



r 



Similarly to the linearizing approximation, this 7^ at the 
UVFP oi ga is consistent with the renormalization oiipip. 
After the renormalization, the effective potential for 
has no divergent term. For < 



,PE 



< Kd < k' 

.PE 



D 



Kd < I'^D^ '-'^ ^^'^ other hand, the four-fermion (crjj) and 
eight-fermion (cr|j) terms diverge in the continuum limit. 



Therefore the corresponding counter terms are required. 
In passing, the bifurcation solution leads to the same 
results as above ones for k^ < kd < k^'* ™ ^^^ whole 
region. Thus only the counter term of the auxiliary field 
propagator is needed in the nonrcnormalizable region < 
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